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ABSTRACT 


The  sky  survey  made  by  the  infrared  astronomical  satellite  (IRAS)  in  1983 
included  observations  of  artificial  earth  satellites.  The  data  base  (with  celestial  objects 
removed)  was  correlated  with  the  NORAD  space  catalogue  to  identify  452  satellites  in 
orbit  above  the  IRAS  900-km  altitude.  The  flux  density  in  three  of  the  four  wavelength 
bands  has  been  analyzed  to  determine  the  temperature,  emissivity,  and  absorptivity  of  the 
identified  resident  space  objects. 
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1.  INTRODUCTION 


Optical  observations  of  artificial  earth  satellites  have  been  made  since  the  launch  of  Sputnik.  Early 
measurements  made  using  reflected  visible  light  [1]  were  for  information  about  the  satellite  position. 
Later,  photometric  observations  were  made  to  investigate  properties  of  the  satellite.  Initially,  signature 
data  were  obtained  in  the  form  of  a  time  series  or  light  curve.  Later,  calibrated  photometric  observations 
were  used  to  discuss  the  reflective  properties,  such  as  reflectivity,  and  configuration  of  satellites  [2]. 
Observations  of  the  self-emitted  IR  radiation  from  satellites  were  also  obtained.  IR  data  are 
complementary  to  the  visible  band  and  UV  data  and  can  provide  information  about  the  absorptivity  and 
emissivity  of  satellite  materials.  These  ground-based  observations  were  frustrated  by  large  and  variable 
atmospheric  absorption  [3]  in  addition  to  sensor  calibration  problems.  Recently,  high-quality  ground- 
based  IR  satellite  observations  have  been  made  [3]. 

The  infrared  astronomical  satellite  (IRAS)  was  launched  in  1983  to  perform  an  all-sky  survey  of  the 
IR  portion  of  the  spectrum.  While  stellar  objects — stars,  nebulae,  comets,  and  asteroids — were  the 
primary  objects  of  interest,  IRAS  also  observed  satellites  and  space  debris.  A  number  of  attempts  have 
been  made  to  extract  these  data  and  demonstrate  that  IRAS  data  can  be  used  to  characterize  IR  emission 
of  orbital  debris.  De  Jong  and  Wesselius  [4]  and  Anz-Meador  et  al.  [5]  used  unprocessed  IRAS  detector 
data  to  search  for  debris.  These  studies  found  a  few  moving  sources  from  searching  only  a  few  days  of  the 
mission  data.  Dow  [6,29]  used  objects  found  in  the  Sky  Brightness  Images  over  the  entire  IRAS  mission. 
All  studies  concluded  that  IRAS  had  the  capability  to  observe  quite  small  debris — Anz-Meador  claimed 
sizes  down  to  1  mm — but  none  systematically  attempted  to  isolate  data  on  known  satellites. 

The  Astronomical  Group  at  Groningen  reprocessed  all  the  raw  IRAS  data  tapes  with  the  objective 
of  finding  all  satellite  and  debris  detections  [4,7,8].  This  effort  produced  a  debris  data  base  containing 
more  than  190,000  detections.  As  a  first  step  in  analyzing  this  data  base,  we  have  concentrated  on  the 
catalogue  of  known  satellites  in  1983.  By  correlating  positions  in  the  debris  data  base  with  the  known 
satellite  catalogue,  we  have  extracted  the  true  satellite  detections  and  been  able  to  discuss  the  associated 
radiometric  data.  This  analysis  is  the  subject  of  this  report. 

The  IRAS  satellite  made  many  observations  of  artificial  earth  satellites.  Processing  for  astronomical 
objects  (stars,  planets,  comets,  and  asteroids),  Beichman  [9]  eliminated  all  satellite  and  debris  data.  The 
Astronomical  Group  at  Groningen  reprocessed  all  the  raw  IRAS  data  tapes.  They  abandoned  the  original 
approach  of  hours,  days,  and  weeks  confirmation,  as  this  method  was  intended  to  eliminate  observations 
of  objects  that  did  not  reappear  at  the  same  place  in  the  sky — just  those  objects  of  interest  here.  They 
established  a  method  that  examined  each  detection  for  multiple  detector  “hits”  that  were  consistent  in 
timing  of  their  passage  across  the  focal  plane.  The  approach  was  very  permissive,  having  wide  bounds  for 
accepting  a  plausible  detection.  The  final  catalogue  contained  more  than  190,000  detections.  The  flux 
values  were  limited  to  those  greater  than  0.1  Jy  (W/m^/Hz).  Also,  they  did  not  consider  data  from  the 
lOO-fj.  band.  This  data  base  was  provided  to  Lincoln  Laboratory.  Consideration  was  limited  to  those 
detections  with  a  signal-to-noise  ratio  (SNR)  greater  than  about  three:  a  data  base  of  about  136,000 
detections. 


2.  IRAS  INSTRUMENT 


The  IRAS  operated  from  March  1983  until  cryogen  depletion  in  November  1983.  IRAS  was  in  a 
sun-synchronous  orbit  with  the  orbital  plane  (800-km  altitude,  99.2°  inclination)  nearly  perpendicular  to 
the  sun  line,  as  illustrated  in  Figure  1 .  The  telescope  was  always  pointed  away  from  the  earth  center,  i.e., 
toward  the  zenith.  The  spacecraft  attitude  control  system  allowed  pole-to-pole  scans  in  ecliptic  longitude 
from  60°  to  1 20°;  most  scans  were  taken  between  84°  and  96°. 
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Figure  1.  IRAS  orbit  geometry. 


The  IRAS  detectors  covered  four  long  wavelength  IR  bands  having  flux  density  reporting  values  of 
12,  25,  60,  and  100  These  correspond  to  peak  flux  densities  of  242,  1 16,  48,  and  29  K.  The  peak 
wavelength  for  the  nominal  free-space  temperature  of  300  K,  7  =  2898/ A,  is  9.66  /i  near  the  lower  cutoff 
of  the  12-/iband. 


3 


WAVELENGTH  BANDS,  pm 


249378-2 


100  60  25  12  60  25  12  100 


\ 

-V . 

I  \ 

I  \ 

■/  ' 

I^CPC  \ 

\ 

I 

I 

/ 

/ 

i  / 

;  / 

■  */ - 

/ 

/ 


LRS 


IMAGE 

DIRECTION 

- 


IRAS  FOCAL  PLANE 


Figure  2.  IRAS  focal  plane. 


The  IRAS  detectors  in  each  of  the  bands  were  arranged  on  a  30-arcmin  focal  plane,  as  illustrated  in 
Figure  2.  A  star  or  satellite  crossing  the  scanning  focal  plane  would  traverse  a  number  of  detectors  in 
several  wavebands.  Each  of  the  59  detectors  had  a  different  response  function,  as  tabulated  in  Beichmann, 
Table  II.C.5,  page  11-18.  The  detectors  had  a  rectangular  shape:  0.75  x  4.5  arcmin  for  the  12-  and  25- fi 
detectors  and  3x5  arcmin  for  the  60-  and  lOO-Ji  detectors.  As  shown  in  Figure  2,  the  spacecraft  attitude 
was  such  that  in  the  scan  direction,  motion  of  fixed  stellar  sources  or  satellites  across  the  focal  plane  was 
along  the  short  detector  axis.  Consequently,  the  positional  accuracy  in  the  along-scan  direction  was  much 
better. 
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The  camera  electronics  used  a  6-Hz  lowpass  filter  to  “despike’’  the  data.  The  detectors  were 
sampled  at  16  Hz.  The  scanning  speed,  determined  by  the  IRAS  orbital  motion,  was  3.85  arcmin/sec. 
Therefore,  a  stationary  source  traversed  the  0.75-arcmin  detector  in  0.19  sec,  giving  about  three  samples 
for  a  detection.  Consequently,  all  objects  (stellar  sources,  solar  system  objects,  satellites,  and  debris)  have 
very  similar  detections.  This  is  illustrated  in  Figure  3  (from  Wesselius  et  al.).  Therefore,  the  response 
function  is  not  a  suitable  vehicle  for  identifying  satellites  and  debris. 
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Figure  3.  Response  due  to  celestial  objects. 


The  baseline  IRAS  method  of  seconds,  hours,  days,  and  weeks  confirmation  for  object 
identification  of  a  detection  relies  on  the  object  being  stationary.  Since  rapidly  moving  satellites  and 
debris  would  be  screened  out,  this  method  was  abandoned.  Wesselius  et  al.  used  the  following  approach. 
Satellites  at  a  geometric  range  greater  than  10,000  km  will  have  an  absolute  velocity  between 
0.25  arcmin/sec  (geosynchronous)  and  10  arcmin/sec  and  apparent  velocity  between  4  and  10  arcmin/sec. 
The  object  can  take  more  than  10  sec  to  cross  the  IRAS  focal  plane  and  will  cross  each  detector  at  a 
different  time,  as  illustrated  in  Figure  4.  Here,  the  detector  samples  are  displayed,  as  a  function  of  time, 
for  detectors  arranged  as  they  were  on  the  focal  plane.  By  correlating  the  time  history  of  detections,  the 
RSO  and  debris  hits  were  selected  from  the  raw  data.  From  this  data  one  obtains  the  amplitude  of  the  hit 
in  each  detector,  the  position,  and  relative  velocity  of  the  object.  Because  of  the  relatively  inaccurate 
position,  the  relative  velocity  is  particularly  important  in  correlating  the  hits  with  a  catalogue  of  known 
satellite  orbits. 
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3.  ANALYSIS  CONSIDERATIONS 


The  IRAS  radiometric  data  is  provided  in  Janskys  uncorrected  for  color.  The  flux  in-band  was 
converted  to  Janskys  “without  prejudice”  by  assuming  a  A”*  frequency  distribution  of  the  objects’ 
radiation.  This  allowed  processing  the  data  without  making  assumptions  about  the  temperature  or 
wavelength  distribution  of  the  source,  and  the  bands  could  be  processing  independently.  Data  reduced  in 
this  way  is  said  to  be  uncorrected  for  color.  The  IRAS  Explanatory  Supplement,  page  11-27,  describes  how 
to  make  the  color  correction  and  obtain  “true”  observed  Janskys.  That  procedure  is  implemented  here 
using  the  tables  given  on  page  VI-26.  Note  that  these  tables  can,  in  principle,  be  derived  from  the  spectfal 
response  function,  Table  II.C.5,  page  11-18.  These  tables  have  been  reconstructed  to  two  significant 
figures.  The  source  of  disagreement  is  unknown,  but-  it  probably  derives  from  subtleties  in  the 
implementation  of  the  calibration  process.  The  basic  process  for  making  the  color  correction  is  as  follows. 

The  observed  flux  density,  at  a  reference  wavelength  fora  source  at  temperature  7,  is 


F  o  (T) 
obs 


,T) 


nr 


where  the  Planck  flux  density  is 


3.74185  xlO’ 


14388.3 


AT 


-1 


w 


\  ^ 


and  X  is  in  microns  (Allen),  A  is  the  projected  area  of  the  object,  e  is  the  emissivity,  and  r  is  the  range  to 
the  object.  The  IRAS  data  is  reported  in  Janskys,  viz.,  10"^^  W/mVHz.  Using  the  relation  that  fX  - 
2.99792458  x  10^"^  ///sec,  one  can  convert  to  Janskys  with 


j\  {  T)  =F^°  {  T) 

ohs  '  '  obs 


1 

'  .  W  ^ 

2.99792458  xlO'^^  j 

^  Hz j 

The  temperature  of  an  object  can  be  found  as  follows.  The  temperature  dependence  of  the  observed 
flux  depends  on  A^.  The  observed  flux  also  depends  on  e.  A,  and  r.  Assuming  e  is  independent  of 

A,  the  ratio  of  observations  made  at  the  same  time 


j\  (T) 

gt(T)=- 


for  bands  A,  and  is  a  monotonic  function  of  T.  Given  a  pair  of  true  flux  densities,  the  temperature  is 
obtained  from  this  function.  Now  we  define  the  color  correction,  Kx,  to  convert  the  quoted,  i.e.,  reported, 
flux  density  fx  =  to  the  “true”  observed  flux  density  as 
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or  converting  from  Janskys  to  W/(m^/x)  then 


F^°  =F^VK 

obs  Q  X 


The  ratio 

9t„(T)=— ^ 

Q 

and  the  Kxs  are  also  monotonic  functions  of  T,  They  can  be  calculated  from  the  spectral  response  function 
for  each  band.  This  calculation  is  described  in  the  IRAS  Explanatory  Supplement,  page  VI-27.  The 
temperature  from  the  flux  density  ratios  /, 2/725  /25//50  can  be  obtained.  The  temperature  obtained  from 
the  quoted  flux  density  ratios  is  the  same  as  that  obtained  from  the  color-corrected  “true”  observed  flux 
densities. 

Initial  modeling  of  RSO  flux  and  temperature  for  experiment  planning  is  done  assuming  that  the 
solar  absorptivity  a  is  equal  to  the  LWIR  emissivity  €.  While  this  is  appropriate  for  solar-cell-powered 
payloads  of  principal  interest,  many  objects,  particularly  debris,  may  have  surfaces  with  the  value  of 
absorptivity  different  from  the  emissivity,  resulting  in  a  large  range  of  free-space  temperatures. 

To  aid  in  the  discussion  a  simple  model  for  the  radiant  flux  is  developed.  This  calculation  requires  a 
number  of  assumptions.  First,  the  object  is  in  thermal  equilibrium,  and  all  the  absorbed  solar  radiation 
energy  is  reemitted  from  each  element  following  Lambert’s  law.  Consequently,  we  ignore  the  energy  that 
is  reradiated  as  microwave  radiation.  Second,  we  assume  there  is  no  contribution  for  earth-upwelling 
radiation.  Third,  we  assume  that  the  satellite  has  one  surface  material,  which  absorbs  solar  radiation  and 
radiates  at  a  single  temperature.  Finally,  we  assume  that  the  satellite  is  not  in  or  near  eclipse.  We  have 
developed  four  simple  models,  which  are  listed  in  Table  1.  These  constitute  four  idealized  satellite 
configurations. 


TABLE  1 

Free-Space  Temperature  for  Basic  Satellite  Shapes 


Principal  Shape 

Temperature 

n=  A,/Ap 

Sphere 

278  K 

4 

Cylinder  (Radiusj;,  -Length  y 

295  K 

Tt(1+  r/L) 

Flat  Solar  Panel  (Two  Sides) 

330  K 

2 

Flat  Solar  Panel  (One  Side) 

393  K 

1 
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The  free-space  equilibrium  temperature  is  determined  by  setting  the  solar  power  absorbed  by  the 
projected  area  equal  to  the  total  power  reradiated  by  the  surface  area  according  to  the  Stephan-Boltzmann 
expression,  in  MKS  units; 

1360  a  A  =  5. 6692x10'®  e  A  T*  . 

P  S 


where  is  the  projected  area  and  >4^  is  the  surface  area.  *  This  relation  then  gives 


T=(393  K) 


For  these  simple  models  using  the  solar  radiation  input  of  1360  W/m^,  one  can  calculate  the  free- 
space  temperature  Ta  =  e,  assuming  that  the  absorptivity  a  equals  the  emissivity  e.  If  the  absorptivity  is 
greater/lesser  than  the  emissivity,  the  temperature  will  be  greater/lesser  than  this  calculated  temperature. 
From  the  observed  temperature  and  emissivity,  we  can  obtain  the  absorptivity  from 


/  ^  4 

a=e  - 

T 

V  a  = 

The  expression  for  total  radiant  intensity  is  then 

1360  a  A 

watts  / steradian  =  - - 

n  n 


where  n  =  A/A^,  as  given  in  Table  1.  Thus,  the  total  radiant  intensity  for  a  sphere  is  lOSoA^, 
watts/steradian.  This  simplified  model  requires  that  the  object’s  radiance  vs.  wavelength  be  a  Planck 
function  with  LWIR  emissivity  e  independent  of  wavelength.  For  objects  in  near-earth  orbit,  the  effects  of 
earth  temperature,  earth-reflected  sunlight,  and  earthshadowing  must  be  included  in  the  RSO  temperature 
estimate  [12]. 


i 


2 

For  this  analysis,  the  value  of  the  solar  constant,  1360  W/m  ,  is  taken  from  Allan,  page  169  [10]. 

2 

measurements  from  Nimbus-7  suggest  a  value  of  1372  W/m  for  1983  [13]. 


Recent 
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4.  SPACECRAFT  MATERIALS 


All  model  calculations  depend  on  the  physical  properties  of  the  satellite  surface  materials.  For  high- 
fidelity  modeling,  detailed  information  is  necessary.  Because  this  information  is  unavailable,  illustrative 
values  are  used,  as  compiled  by  Dow  [6]  and  reproduced  in  Table  2  for  a  and  e.  Also  given  are 
temperatures  for  a  sphere^  (Tl)  and  flat  plate  (T2)  covered  with  the  material.  The  values  of  ale  range 
from  0.1  to  10.0.  Figure  5  plots  reflectivity  (p  =  1  -  a)  vs.  the  tabulated  ale  ratios.  As  seen,  almost  all 
values  of  a  and  e  are  possible.  For  example,  the  temperature  of  a  gold-plated  sphere  can  rise  above  the 
nominal  278  K,  for  a/6  =  1,  to  494  K,  and  a  sphere  plated  with  magnesium  oxide  white  paint  can  drop  to 
156  K.  Single-sided  flats  can  rise  from  the  nominal  393  K,  for  ale  =  1,  to  about  700  K! 


Figure  5.  Reflectivity  vs.  o/e  for  spacecraft  materials. 


There  is  a  paucity  of  data  on  the  dependence  of  e  with  variables  such  as  wavelength  and 
temperature.  Ross  et  al.  [13]  report  measurements  on  silicon  solar  cells  where  they  give  the  value  of  a  = 
0.78,  and  emissivity  changes  from  0.68  at  -50®C  to  0.78  at  1 50®C  (see  Figure  6).  Blair  et  al.  [14]  report  on 
extensive  measurements  on  a  number  of  materials  with  emphasis  on  solar  cells  manufactured  by 
AEROJET.  They  measure  the  emissivity  dependence  on  reflectance  angle,  temperature,  and  wavelength 
ranging  from  2  to  24  //  (see  Figure  7).  From  2  to  18  p  the  emissivity  is  generally  between  0.8  and  0.9. 
There  is  a  decreasing  trend  in  emissivity  for  wavelengths  greater  18  p,  although  the  data  are  inconclusive. 


^  Dow  inverted  this  table.  Dow  used  1400  W/m^  for  the  solar  power,  whereas  1360  W/m^  is  used  in  this 
analysis. 
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They  report  emissivity  of  about  0.9  at  200  K  and  less  at  373  K.  This  is  the  opposite  trend  from  that 
reported  by  Ross  et  al.  Dow  reports  that  the  emissivity  of  white  paint  changes  from  0.9  at  50  fi  to  0.2  at 
75  fx  and  to  0.1  at  100  fi.  Therefore,  some  dependence  of  emissivity  on  temperature  is  expected,  although 
even  the  sign  of  the  slope  is  unknown.  More  probably  we  would  expect  the  emissivity  to  decrease  with 
increasing  wavelength.  This  is  consistent  with  Drudes’s  theoretical  relation  between  the  resistivity  q  of  a 
metal  for  a  definite  wavelength  [15]. 
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Figure  6,  Emissivity  vj.  temperature  for  silicon  solar  cells. 
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Figure  1.  Emissivity  vs.  wavelength  for  silicon  solar  cells. 
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Currently,  it  is  not  generally  known  if  the  values  of  ale  change  with  time  in  space.  For  optical  solar 
radiators,  ale  does  change  from  0.08  to  0.2  or  more  over  10  years  [18].  This  is  inferred  from  the  increase 
in  electronic  component  temperature.  Results  from  the  Long  Duration  Experiment  Facility  (LDEF),  in 
orbit  for  almost  six  years,  provide  some  indication.  As  summarized  by  Dow,  for  anodized  aluminum  (the 
primary  LDEF  surface  coating),  the  values  of  e  did  not  change  significantly,  whereas  the  value  of  a 
increased  by  as  much  as  16%,  depending  on  the  location  of  the  sample.  However,  white  paint  exhibited 
increases  as  large  as  1 00%.  Such  phenomenon  may  help  explain  the  high  and  low  temperatures  reported 
here. 
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TABLE  2 

Typical  Absorptivity  and  Emissivity  for  Spacecraft  Materials 


Material 

P 

a  =  1-p 

e 

a/e 

T1 

T2 

AlO,  on  Buffed  Aluminum 

0.87 

0.13 

0.23 

0.57 

243 

289 

Vapor  Dep.  Gold  on  Glass 

0.96 

0.04 

0.02 

2.00 

333 

396 

Vapor  Dep.  Silver  on  Glass 

0.81 

0.19 

0.02 

9.50 

492 

585 

Solar  Cell  (lUE) 

0.14  ■ 

0.86 

0.84 

1.02 

282 

335 

Stainless  Steel  (Polished) 

0.58 

0.42 

0.11 

3.82 

392 

465 

Fiberglass 

0.15 

0.85 

0.75 

1.13 

289 

344 

Polished  Aluminum  6061 

0.81 

0.19 

0.04 

4.52 

408 

485 

Unpolished  Aluminum  6061 

0.63 

0.37 

0.04 

8.81 

483 

574 

Anodized  Aluminum 

0.58 

0.42 

0.63 

0.67 

254 

302 

Silver 

0.96 

0.04 

0.02 

2.00 

333 

396 

Gold  (Plated) 

0.70 

0.30 

0.03 

10.0 

498 

593 

TiO  White  Paint 

0.80 

0.20 

0.90 

0.22 

192 

228 

Black  Paint(#M  Velvet) 

0.05 

0.95 

0.92 

1.03 

282 

336 

Aluminum  Paint 

0.71 

0.29 

0.27 

1.07 

297 

353 

1 .0  mil  Aluminized  Mylar* 

0.84 

0.16 

0.54 

0.30 

207 

247 

1.0  mil  Silverized  Teflonf 

0.92 

0.08 

0.66 

0.12 

165 

196 

1 .0  mil  Aluminized  Kaptont 

0.64 

0.36 

0.54 

0.67 

254 

302 

Magnesium  Oxide  White  Paint 

0.91 

0.09 

0.90 

0.10 

158 

187 

Platinum  Foil 

0.67 

0.33 

0.04 

8.25 

475 

565 

*  Trade  name  for  polyethylene  terephthalate 

t  Trade  name  for  fluorinated  ethylene  propylene 

t  Trade  name  for  polyamide 
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5.  CALIBRATION  ISSUES 


The  IRAS  data  received  from  Groningen  was  processed  with  the  calibration  described  in 
Beichman’s  report  [9].  Since  that  publication,  a  number  of  analyses  have  questioned  the  calibration 
[16,17,18,19,20,21].  Some  consider  the  calibration  confirmed  (Aumann  [22]),  whereas  a  number  of 
analyses  suggest  the  various  wavebands  are  in  error  by  2%  to  12%.  The  proposed  corrections  are  listed  in 
Table  3.  The  weight  of  evidence  suggests  that  some  calibration  correction  may  be  warranted.  The 
calibration  correction  adopted  for  this  analysis  is  also  given  in  Table  3.  In  this  analysis,  there  is  no  way  to 
quantify  or  assess  the  absolute  calibration,  as  we  do  not  have  independent  information  about  the  satellite 
emissivity,  absorptivity,  or  temperature.  It  is  the  determination  of  these  quantities  that  is  the  objective  of 
this  analysis. 


TABLE  3 

IRAS  Calibration  Errors 


Source 

12;/ 

25 

60  ^J 

100  tJ 

Aumann 

0 

0 

0 

0 

Tedesco 

0 

10% 

10% 

7 

Gillett 

0 

6% 

0 

7 

Kirby  et  al. 

0 

10% 

0 

0 

Cohen  et  al. 

2% 

6% 

3% 

12% 

Used  Here 

0 

0 

0 

NA 
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6.  DATA  PROCESSING 


Each  detection  report  contains,  among  other  information,  the  epoch,  direction  in  right  ascension  and 
declination,  angular  rate  crossing  the  focal  plane,  and  the  flux,  in  uncolor-corrected  Janskys,  for  each 
detector  registering  a  “hit.” 

The  objectives  of  this  analysis  are  to  identify  those  detections  of  known  satellites  and  to  analyze  the 
radiometric  data  for  information  about  IR  observations  of  artificial  earth  satellites.  There  were  four  steps 
in  preparing  the  data  for  analysis. 

1.  The  observed  direction  of  all  detections  in  the  Groningen  debris  data  base  was  compared  with 
directions  computed  from  an  ephemeris  of  the  IRAS  satellite  and  the  catalogue  of  known  satellite  element 
sets.  Observations  were  accepted  if  the  observed  and  computed  direction  agreed  to  0.60°  or  better  and  the 
observed  and  computed  angular  velocity  agreed  to  0.5  arcmin/sec. 

2.  The  individual  hits  for  each  detection  were  screened.  The  intention  was  to  eliminate  hits  where 
the  satellite  image  (assumed  a  point  source)  did  not  completely  cross  a  detector,  as  a  partial  hit  would  give 
a  biased  flux  estimate.  The  screening  was  done  as  follows.  As  seen  from  the  focal  plane  geometry  (see 
Figure  2),  the  detectors  for  each  waveband  are  in  two  rows  with  some  overlap.  Therefore,  in  each 
detection  one  would  expect  that  two  hits  would  be  registered  for  each  waveband.  In  cases  where  three  or 
more  hits  are  registered,  it  is  assumed  that  some  of  them  are  at  a  detector  edge  and  therefore  do  not  give  a 
good  flux  measurement.  Therefore,  the  screening  was  that  for  detections  with  three  or  more  hits,  the  two 
hits  with  the  largest  flux  density  measurement  were  accepted.  This  eliminated  some,  but  probably  not  all, 
partial  detector  hits.  An  example  of  the  editing  of  the  data  for  SBSl  (12065)  is  presented  in  Table  4.  The 
sighting  time  for  each  detection  is  given  as  the  day  of  the  year  1983.  The  record  number  is  a  unique 
identifier  in  the  Groningen  IRAS  debris  data  base.  The  object  number  is  12065.  The  observed  in-scan 
velocity  is  given  in  arcmin/sec.  In  each  detection,  all  the  hits  are  given,  and  the  detector  is  identified.  A 
number  of  triple  hits  are  observed.  The  discarded  hit  is  noted  with  an  asterisk.  Hits  were  discarded  for 
other  reasons.  The  detection  on  day  222.032  was  discarded  because  of  the  anomalous  in-scan  velocity. 
Other  hits  discarded,  based  on  a  statistical  test,  were  in  detections  on  days  71.535  and  72.537. 

3.  The  detections  for  each  satellite  were  collated. 

4.  The  flux  density  measurements  were  color  corrected  and  analyzed  to  obtain  physical  properties 
of  the  satellite  using  other  observed  information  such  as  the  range  and  physical  size.  Information  such  as 
temperature  and  emissivity  were  derived. 
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TABLE  4 

Detector  Hits  for  SBS1 
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7.  GENERAL  RESULTS 


Table  5  summarizes  the  result  of  correlating  with  the  satellite  catalogue.  The  number  of  correlations 
is  given  for  a  variety  of  error  bounds.  The  2072  detections  with  angle  errors  less  than  0.6°  are  believed  to 
be  valid  satellite  detections.  The  remaining  detections  contain  uncatalogued  objects  (UCTs),  debris,  and 
false  detections.  These  will  not  be  discussed  further. 


TABLE  5 

IRAS  Data  Correlation  Results 


Number  of  IRAS  sightings  compared  with  catalogue 

136304 

Number  of  correlations 

Crude  filtering 

30600 

<4‘ 

9955 

<  r 

2469 

<0.6' 

2072 

Number  of  satellites  correlated 

465 

Statistics 

vs.  IRAS  DB 

vs.  Correlations 

Correlations 

22% 

<4* 

7% 

33% 

<r 

2% 

8% 

<0.6’ 

2% 

7% 

The  Wesselius’s  original  publication  [8]  displayed  the  data  in  an  interesting  way.  From  the 
observed  flux  ratio,  say  one  can  determine  the  temperature  of  the  body,  assuming  the  emissivity 

is  constant.  One  would  expect  the  temperature  measured  with  the/|2;i//25p  ratio  to  be  the  same  as  that 
obtained  from  ratio.  This  is  plotted  in  Figure  8  from  Wesselius  et  al.  for  the  1500  detections 

with  fluxes  greater  than  3  Jy.  The  wide  scatter  belies  our  expectation.  The  same  plot  is  provided  for  the 
correlated  data  before  the  screening  for  partial  hits  (see  Figure  9).  Instead  of  selecting  detections  based  on 
large  flux  density,  we  have  selected  detections  based  on  known  satellite  detections.  A  similar  scatter  is 
evident.  Next,  the  same  plot  is  provided  for  temperature  derived  from  the  screened  flux  density 
measurements  (see  Figure  10).  In  this  case,  the  values  fall  on  a  straight  line,  and  the  scatter  is  significantly 
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reduced.  From  the  f\nJf2iti  flux  ratio,  the  average  temperature  of  about  300  K  is  clearly  evident.  We 
believe  the  large  scatter  in  the  original  IRAS  data  is  due  in  part  to  partial  detections.  In  principle  a  partial 
detection  could  be  corrected  given  precise  knowledge  of  the  IRAS  position  and  pointing  and  the  satellite 
position. 

Another  way  to  characterize  the  IR  satellite  measurements  is  given  in  Figure  1 1.  Here,  the  12-/rflux 
density  is  presented  as  a  function  of  the  range  to  the  observed  target.  The  envelope  corresponds  to  the 
dependence,  and  its  position  depends  on  the  sensitivity  of  IRAS  and  the  maximum  size  of  objects  in  1983. 
A  number  of  detections  are  at  a  greater-than-geosynchronous  range.  These  are  genuine  since  they  are 
Astron,  Exosat,  and  other  known  satellites.  The  individual  fluxes  are  given  in  Appendix  A. 
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Figure  8.  Temperature  from  1500  large  flux  sample  detections. 
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Figure  9.  Temperature  from  correlated  and  unscreened  data. 


Figure  10.  Temperature  from  correlated  and  screened  data. 
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Figure  1 1.  Flux  density  vs.  range. 
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8.  REFERENCE  SPHERES 


Isothermal  spheres  are  the  simplest  and  easiest  satellite  shapes  to  analyze.  A  single  IRAS  focal 
plane  crossing  has  been  correlated  with  the  polished  aluminum  Lincoln  Radar  Calibration  Sphere  LCSl 
(SSC  #  1361).  Although  there  is  only  one  detection,  the  analysis  of  this  data  is  instructive.  The 
temperature  is  determined  to  be  480  K  using  the  screened  uncolor-corrected  flux  density  ratio  of  = 

17.69  Jy  at  12  /t  over  1 1  22  Jy  at  25  /t.  The  radiant  flux  density  of  LCSl  at  480  K  is  = 

134.8  W//t/(effective  eA)  at  a  12-/t  reference  wavelength. 

The  IRAS  color  correction  (Beichmann,  Table  VI-26)  adjusts  the  reported  M-fi  flux  density  to 
=  1 6.55  Jy,  which  is  equal  to  =  F^'^^‘IK^2f^  =  3.45  x  10’'^  W/(m^/t)  spectral  irradiance.  Now  at 

a  range  of  2001  km,  determined  from  the  a  posteriori  ephemeris,  the  observed  spectral  radiant  intensity  is 

;i^F,^i2#i  =  4.34  w//i. 

Now  the  effective  area  eA  is  the  observed  spectral  radiance  at  12  (W//t)  divided  by  the  Planck 

function  value  of  the  radiant  flux  density  [W/(m^  fi)]  at  the  same  reference  wavelength. 

nr2pi2u 

pA  =  _ 


0.0322  (m^) 


Thus,  eA  is  0.0322  m^,  but  since  the  projected  area  for  LCSl  (1361),  A  =  1  m^,  the  emissivity  is  e  = 
0.0322. 

The  ratio  of  solar  absorptivity  to  emissivity  (a/e)'^  is  equal  to  the  observed  temperature,  480  K, 
divided  by  278  K,  the  characteristic  free-space  temperature  for  an  isothermal  sphere.  Thus,  a/e  =  8.89, 
a  =  0.286,  and  the  reflectivity  is  p  =  1  -  a  =  0.7 1 39. 

Data  on  the  optical  properties  of  various  spacecraft  materials  give  values  for  aluminum  as  shown  in 
Table  6. 

The  radar  calibration  sphere  may  have  a  less  polished  surface  than  when  launched  a  quarter  century 
ago,  thereby  increasing  the  absorptivity — i.e.,  reducing  the  reflectivity — toward  the  values  given  for 
unpolished  aluminum. 

The  IRAS  result  can  be  compared  with  visible  band  measurements.  The  apparent  visual  brightness 
resulting  from  a  specular  sphere  with  pA  =  0.714  m^  at  a  2(X)0-km  range  is  given  by 


M =-26.78-2 . 5 


/ 

log 


0.714  ^ 

4n  (2x10^)^  j 


7 . 84M 

V 
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TABLE  6 

Optical  Properties  of  Aluminum 


a 

e 

ale 

P 

T(K) 

LCSl  Observed 

0.286 

0.0322 

8.89 

0.714 

480 

Polished  Al 

0.19 

0.042 

4.52 

0.81 

408 

Unpolished  Al 

0.63 

0.042 

8.81 

0.63 

483 

This  calculated  brightness  can  be  compared  with  the  brightness  observed  by  W.  Beavers  [2]  as  a  function 
of  phase  angle.  The  value  derived  from  IRAS  for  is  7.84  at  2001  km.  This  can  be  referred  to 
geosynchronous  distance,  =  36,000  km,  with  an  adjustment  of  6.275  Af,,.  Therefore,  IRAS  predicts  14.1 1 
at  geosynchronous  range.  This  is  in  moderate  agreement  with  Beavers’s  observations  reported  at 
geosynchronous  range  that  give  a  value,  at  90°  phase  angle,  of  14.0  M^.  Beavers  notes  that  his 
observations  indicate  that  12%  of  the  energy  is  diffuse.  Additional  LWIR  and  visible  band  measurements 
at  the  same  epoch  are  required  for  a  high-accuracy  comparison. 

Finally,  consider  the  thermal  response  of  LCSl,  a  1-m  sphere,  with  a  2-cm-diameter  solid 
aluminum  sphere  covered  with  Martin  Black  with  a  =  e  =  0.999  as  it  enters  the  earthshadow.  The  cooling 
rate,  neglecting  earth-upwelling  radiation,  is 

dr__  5.67x10'^  e  A  (1-^)  [  Jcw' 

dt  c^M  I  ’ 

where  is  the  specific  heat,  M  is  the  mass,  and  (1  -  is  the  fraction  of  the  sphere  radiating  to  free  space. 
Now,  the  temperature  of  a  sphere  is  written  as 


T= 


278 


\ 
a 


giving 


dr__  338.65  g  A  (l-Q 
dt  c  M 

p 
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Now,  the  area-to-mass  ratio  of  the  two  spheres,  (A/M)(,35|)  =  0.0285  m^/Kg,  and  (A/M)(2  cm)  = 
0.0252  m^/Kg,  is  virtually  the  same.  The  cooling  rate  for  LCSl  (1361)  is  therefore  expected  to  be  about 
0.286  that  of  the  2-cm  emissive  calibration  sphere  coated  with  Martin  Black,  which  results  in  a 
temperature  excursion  less  than  10  K.  Therefore,  LCSl  would  provide  a  relative  stable  temperature 
reference. 
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9.  SPIN-STABILIZED  CYLINDERS 


Spin-stabilized,  solar-cell-covered  cylinders  appear  to  be  simple  objects  for  radiometric  analysis. 
There  were  80  such  nonmilitary  satellites  in  orbit  during  the  IRAS  sky  survey  compared  to  20  three-axis 
stabilized  satellites.  The  putative  simplicity  of  spin-stabilized  cylinders  led  to  early  analysis  of  that  data. 

The  first  geostationary  satellite,  Syncom  B,  a  28-in-diameter  and  15-in-high  spun  cylinder,  was 
launched  in  1963.  The  cylinder  weighed  85  lb.  The  solar  power  was  30  W,  and  the  satellite  had  a  lifetime 
of  18  months.  Now,  30  years  later,  Hughes  (the  same  company  that  built  Syncom)  deploys  two-ton 
cylinders,  1 2  ft  in  diameter  and  30  ft  high  with  2000  W  of  solar  power  at  the  end  of  a  ten-year  life. 

The  first  cylinders  employed  dipole  antennas,  sometimes  in  phase  switched  arrays,  to  provide 
moderate  gain  in  the  earth  direction.  Later,  satellites  used  high-gain  parabolic  reflectors  on  a  despun 
platform  at  the  top  of  the  cylinder.  The  other  end  of  the  cylinder  was  left  open,  and  a  second  concentric 
solar-cell-covered  cylinder  slid  down  over  the  apogee  kick  motor  after  geosynchronous  orbit  was 
achieved.  These  geosynchronous  orbits  have  nearly  0®  inclination  and  circular  orbits  with  the  spin  axis 
maintained  normal  to  the  orbit  plane,  i.e.,  normal  to  the  earth’s  equatorial  plane.  The  IRAS  orbit  geometry 
results  in  observations  of  geosynchronous  satellites  near  quadrature,  i.e.,  at  an  illumination  phase  angle 
between  80*"  and  100®.  The  measured  flux  density  from  these  cylinders  is  exi:)ected  to  be  a  weak  function 
of  phase  angle,  which  is  in  contrast  to  flat  solar  panels  that  exhibit  a  cos(0)  relationship.  The  LWIR 
detections  of  spin-stabilized  cylinders  are  therefore  expected  to  be  the  simplest  to  understand. 

The  simple  LWIR  mode!  used  here  treats  only  the  main  cylindrical  body.  The  LWIR  radiance  of  the 
antenna  and  despun  platform  are  ignored  for  two  reasons.  First,  they  are  covered  with  reflective  thermal 
insulating  material.  Second,  because  of  the  IRAS  geometry,  the  observations  are  made  when  the  parabolic 
antenna  is  edge-on  to  the  sun,  it  is  cool,  and  its  radiation  is  ignorable.  The  LWIR  energy  radiated  from  the 
open  bottom  of  the  cylinder  is  also  ignored  in  the  energy  balance  calculation.  Now,  the  cylinder,  spinning 
on  the  order  of  once  per  second,  is  certainly  near  thermal  equilibrium,  and  only  the  projected  area  is 
considered.  Incidentally,  the  cylindrical  geometry  produces  only  l/;ras  much  power  as  an  equal  number 
of  solar  cells  on  a  flat  pane!  norma!  to  the  sun. 

From  information  about  solar  cell  absorptivity  and  emissivity  it  is  known  that  about  10%-15%  of 
the  incident  solar  power  is  reflected  by  the  solar  cells  p  =  1  -  a.  About  15%  of  the  incident  solar  power  is 
converted  to  electrical  power.  This  leaves  about  75%  of  the  solar  power  to  be  reradiated  to  achieve 
thermal  equilibrium.  The  free-space  temp>erature  of  a  solar  panel  is  changed  by 

where  e  «  0.85,  and  a  =  0.75.  Note  that  a  =  1  -  p  -  7],  where  the  symbol  for  reflectivity  is  p  and  the 
symbol  for  solar  cell  efficiency  is  7].  It  is  the  fourth  root  of  the  ratio  of  this  net  absorptivity,  a  =  0.75,  to 
emissivity,  e  =  0.85,  that  determines  the  temp)erature  of  a  given-shape  solar  panel. 

Most  cylindrical  payloads  were  built  by  Hughes,  including  the  giant  HS-393  class,  which  includes 
SBS-6  and  INTELSAT-6  (all  launched  after  1983).  The  Hughes  Galaxy  class,  HS-376,  Figure  12,  was 
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selected  for  initial  analysis  because  of  the  number  of  IRAS  detections,  «  60,  on  eight  individual  satellites 
that  were  obtained.  In  addition,  they  had  flux  densities  exceeding  2  Jy.  Now,  the  HS-376  class  was  the 
first  to  radiate  heat  from  the  electronics  using  a  cylindrical  optical  solar  radiator  (OSR).  Previous  designs 
radiated  the  heat  from  the  electronics  upward  toward  the  antenna  via  the  despun  platform.  To  obtain  the 
simple  model  for  the  cylinder,  it  is  assumed  that  the  radiator  power  is  equal  to  the  reradiated  power  of  the 
solar  panel  it  replaced.  The  HS-376  was  also  the  first  to  increase  the  solar  cell  collection  area  by  sliding  a 
concentric  solar  cylinder  down  to  cover  the  apogee  kick  motor  after  geostationary  orbit  was  achieved. 
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Figure  12.  Hughes  HS-376  satellite. 
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The  eight  Hughes  cylinders  observed  by  IRAS  are  listed  in  Table  7.  Figure  13(a),  (b),  and  (c)  plots 
the  detections  as  a  function  of  time;  there  is  no  indication  of  a  temporal  trend,  and  we  would  expect  none. 
This  confirms  the  constancy  of  the  IRAS  calibration.  Therefore,  we  have  taken  an  average  of  all  the 
screened  detections  and  obtained  the  mean  values  given  in  Table  7,  which  are  accompanied  by  the  formal 
uncertainty.  In  general  the  average  flux  is  obtained  to  ±  0.1  Jy.  With  this  statistical  population,  the 
accuracy  of  a  single  hit  is  estimated  as  ±  0.6  Jy.  The  object  temperature  is  obtained  as  described  above. 
We  adopted  the  temperature  derived  from  the/,2p//25p  ratio  as  the  object  temperature  for  further  analysis. 
Given  the  temperature,  one  can  then  recover  physical  properties  of  the  satellite  as  follows.  The  observed 
flux  densities  are  color  corrected  to  find  the  true  flux  density.  Since  the  range  is  known  from  the 
ephemeris  used  for  correlation  and  the  self-emitting  projected  surface  area  of  these  satellites,  = 
8.84  m^,  the  emissivity  can  be  calculated  in  each  waveband  from 
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The  calculation  of  emissivity  for  the  eight  Hughes  cylinders  is  given  in  Table  8.  Note  that  the  12- 
and  25- ji  band  emissivities  are  in  good  agreement  with  each  other,  as  they  should  be  since  they  were 
assumed  to  be  so  for  calculation  of  the  temperature.  The  12-  and  25-/i  emissivities  are  also  in  good 
agreement  with  published  values  (Table  2),  whereas  the  60-/i  values  are  not.  We  assume  that  either  there 
is  an  (unlikely)  calibration  error  in  the  60- /i  band  or  there  is  a  real  difference  in  solar  cell  emissivity  at 
60  ji.  Finally,  the  solar  absorptivity  of  the  satellites  can  be  calculated  as  described  above.  The  absorptivity 
calculated  in  this  way  is  also  listed  in  Table  8. 


From  the  HS-376  results,  a  convenient  rule  of  thumb  is  postulated  relating  the  observed  M-ji  flux 
7,2^  to  projected  solar  cell  area 

5  J,.,..  im^) 


12u 


This  leads  to  a  relation  between  power  and  observed  flux 
W  -  500  {W) 

12  p  '  ' 

for  nominal  10%  solar  cell  efficiency. 
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Figure  13.  Flux  vs.  time  for  Hughes  HS-376  cylinders. 
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The  general  result  for  analysis  of  spin-stabilized  cylinders  is  that  relatively  simple  models  can  give 
important  information  about  these  satellites.  Assuming  that  solar  cells  are  the  primary  source  of  LWIR 
flux  seems  a  workable  assumption.  It  allows  inference  of  effective  emissivities  and  possibly  some 
inference  about  the  solar  cell  power  available. 

The  HS-376  spacecraft  observed  by  IRAS  were  assumed  initially  to  be  more  or  less  identical  and 
were  expected  to  yield  the  same  flux  densities  and  temperatures.  The  values  of  flux  density  for  a  given 
spacecraft  were  averaged  without  adjustment  for  possible  seasonal  variations.  However,  individual 
spacecraft  had  important  differences,  some  of  which  could  explain  the  difference  in  observed 
characteristics.  WESTAR  4  and  5  have  shorter  solar  panels  and  corresponding  lower  power  [23].  If  the 
radiator  is  the  same  size,  which  is  not  documented,  the  two  WESTARs  have  only  7.8  m^  projected  area 
compared  with  the  other  HS-376  with  8.84  m^.  Thus,  the  WESTARs  have  0.88  of  the  area,  and  this  is  in 
good  agreement  with  the  1  l-fx  flux  density  ratio  of  0.87  of  the  average  of  the  two  WESTARs  compared 
with  the  two  ANIKS  (1.613  Jy/1.851  Jy).  In  addition,  the  two  SBS  satellites  had  IRAS-derived 
temperatures  about  15  K  higher  than  the  other  spacecraft.  This  may  be  due  to  the  smaller  number  of  20-W 
TWT  amplifiers  employed.  Yenne  [24]  reports  that  SBS  satellites  had  only  10  TWTs  compared  with  the 
usual  20  on  other  satellites.  Another  possible  source  of  difference  is  that  early  spacecraft  were  launched 
without  baking  out  the  solar  cell  panels  to  remove  the  water  from  the  hydroscopic  solar  cell  adhesives. 
Both  SBS-1  and  WESTAR  4  were  not  baked  out  before  launch,  resulting  in  a  few  percent  lower  solar 
power  for  several  months  during  outgassing  in  orbit.  However,  outgassing  should  have  been  completed 
long  before  the  IRAS  data  was  taken. 

The  most  startling  difference  among  spacecraft  is  the  switch  from  a  K7  Spectrolab  solar  cell  to  the 
mixture  of  K7  and  K4%  cells  on  the  Galaxy  1.  The  K7  cells  have  been  featured  as  having  the  highest 
peak  solar  output  (19  mW/cm^),  albeit  with  higher  temperature  resulting  from  a  higher  absorptivity,  a  ~ 
0.84.  The  K7  cells  are  thin  (2.5  mils)  and  have  a  textured  front  surface  together  with  an  aluminized  rear 
coating  that  reflects  near-IR  back  through  the  cell  for  increased  power  generation.  The  K7  was 
particularly  extolled  for  use  on  the  deployable  solar  panel  where  the  characteristic  excess  heat  could  be 
easily  radiated  [25].  Apparently,  the  K7  cells  on  the  fixed  panels  flanking  the  radiator  made  it  more 
difficult  to  control  the  temperature  of  the  electronics  and  batteries  behind  the  radiator  panel. 

The  design  of  high-power,  long-lived,  solar-cell-covered  cylinders  is  strongly  driven  by  the 
temperature  control  required  for  electronic  components  and  batteries,  which  power  the  RSO  during  earth 
eclipse  near  equinoxes.  The  solar  cells  are  cemented  to  1/2-in-thick,  Kevlar-faced,  aluminum  honeycomb 
structural  panels.  In  the  early  Hughes  376  series,  the  new  K7  high-efficiency  solar  cells  were  used  on  both 
the  fixed  upper  (forward)  and  deployable  (lower)  cylinders.  The  switch  to  the  cooler  K4  solar  cells 
suggests  a  temperature  control  issue. 
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The  OSR  has  quartz  panel  sections  with  second  surface  mirrors  and  black  backing  coatings.  The 
initial  solar  absorptivity-to-reflectivity  ratio  is  less  than  0.08.  Thus,  the  radiator  cooling  section,  between 
the  two  upper  solar  panels,  acts  as  a  one-way  mirror  that  reflects  sunlight  and  allows  heat  from  the  TWTs, 
nickel  cadmium  battery  cells,  etc.,  to  radiate  into  free  space.  The  electronics  shelf,  including  the  RF 
output  TWTs  behind  the  radiator,  is  on  the  despun  platform  with  the  antenna  and  feeds.  The  batteries  spin 
with  the  solar  panels.  The  rise  of  spacecraft  temperature  of  several  degrees  over  its  lifetime  has  been 
attributed  to  the  slow  degradation  of  the  OSR’s  a/e. 

The  power  system  for  spacecraft  usually  comprises  two  redundant  main  array  buses  together  with 
battery  charge  and  trickle  arrays.  The  nearly  14,000  solar  cells  of  the  HS-376  spacecraft  are  connected  in 
various  “strings”  of  series-parallel  configurations,  with  protective  diodes,  to  allow  graceful  degradation  in 
the  event  of  string  failure. 

Hughes’s  modeling  of  solar  cell  power,  voltage,  and  current  predicts  performance  within  a  few 
percent.  This  is  confirmed  by  several  years  of  telemetry.  Temperature  is  also  modeled,  but  telemetric  data 
agrees  with  models  to  only  ±6%  of  the  typical  solar  panel  temperature  (283  K).  The  seasonal  variation  of 
bus  voltage  due  to  the  23.5“  obliquity  of  the  ecliptic  and  the  earth  orbit  eccentricity  is  clearly  observed  in 
the  models.  The  temperature  cycle  is  not  so  apparent  since  at  eclipse  season,  when  more  solar  power  is 
available  during  sunlight,  more  power  must  be  used  to  recharge  the  batteries  following  the  eclipse. 

More  detail  about  solar  cell  panels  is  given  for  the  HS-376  successor,  the  HS-393  [26].  Table  9, 
taken  from  Fodor  et  al.  converted  to  degrees  Kelvin,  summarizes  the  predicted  array  panel  temperatures 
by  season  and  for  the  beginning  and  end  of  life.  K4)^  cells  are  used  only  on  the  fixed  panels,  and  the  K7 
cells  are  used  on  the  deployable  panel.  Though  the  total  power  of  the  two  buses  is  over  2000  W,  the 
number  of  solar  cells  required  is  not  larger  than  for  HS-376  because  much  larger  solar  cells  have  been 
developed  to  simplify  fabrication.  Note  that  the  temperature  of  the  cell  panel  is  higher  by  5-9  K 
than  the  K7  panels,  except  near  winter  solstice.  This  indicates  that  equipment  temperatures  are  better 
controlled  with  the  K4  X  cells  near  the  radiator.  Incidentally,  the  lower  panel  falls  to  191  K  before 
emerging  from  eclipse.  Also  note  the  trickle  charge  panel  rises  to  306  K  at  winter  solstice  as  the  sunlight 
partially  illuminates  the  inside  of  the  lower  cylinder. 

It  was  recently  learned  that  in  the  mid-1970s  solar  cell  improvements  foretold  an  operational 
increase  in  watts  per  square  meter  of  30%  to  50%.  Allison  et  al.  [27]  describes  an  improved  “violet”  cell 
(announced  in  1972)  with  vastly  improved  response  to  the  blue-violet  solar  spiectrum  that  achieved  at  least 
a  30%  increase  in  solar  power  over  the  conventional  cell  in  use.  The  COMSAT  nonreflective  (CNR)  cell 
announced  in  1974  offered  a  50%  increase  in  power  above  the  conventional  cell  as  well  as  increased 
radiation  hardness.  The  solar  efficiency  t]  is  15.4%. 
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TABLE  9 

Nominal  Solar  Panel  Temperature  for  HS-393  Satellites 


K4y4 

K7  Deployable 

Date 

Years 

in 

Array  Temperature  (Deg  K) 

Panel 

Panels 

Panel 

Panel 

Orbit 

A1 

C1-C4  MCI 

MC2 

Tc 

Summer 

0 

275.8 

266.9 

268.0 

265.2 

Solstice 

10 

278.0 

266.9 

268.0 

265.2 

Autumnal 

0 

280.8 

275.8 

277.4 

272.5 

Equinox 

10 

283.0 

275.8 

277.4 

272.5 

Winter 

0 

276.3 

297.4 

289.1 

305.8 

Solstice 

10 

278.6 

298 

289.7 

305.8 

21  March 
Post-EcI 

0 

191.3 

144.7 

Meulenberg  [20]  also  advises  that  (a  -  TjVe  is  a  better  indicator  of  cell  temperature  than  ale. 
Table  10,  taken  from  the  COMSAT  paper,  tabulates  the  properties  of  conventional  cells,  the  violet  cell, 
and  the  CNR  cell.  We  have  added  the  calculated  cell  temp)erature  modeled  on  a  cylindrical  satellite 


T  =  295 


a  -  n 


The  CNR  cell  clearly  has  a  higher  temperature,  and  the  paper  suggests  that  a  specific  cover  slide 
would  increase  the  emissivity  from  0.803  to  0.843,  thereby  somewhat  ameliorating  the  temperature 
increase.  The  K7  cells  used  on  the  early  Hughes  HS-376  class  satellites  are  more  or  less  equivalent  to  the 
CNR  cells.  IRAS  data  on  seven  of  the  eight  observed  satellites  were  reexamined — the  PALAPA  B1 
detection  was  excluded — to  correlate  the  observed  temperatures  with  the  solar  cell  models.  The  results  are 
given  in  Table  1 1 . 
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TABLE  10 

Solar  Cell  Types  and  Temperature 


Type 

Electrical 

Efficiency 

Solar 

Absorptivity  a 

Normal 

Emittance 

(a-n)/e 

Temperature 

Conventional 

(Pre-1976) 

0.101 

0.725 

0.803 

0.777 

277 

0.101 

0.725 

0.843* 

0.74 

273.6 

"Violet" 

0.14 

0.81 

0.803 

0.834 

281.9 

0.14 

0.81 

0.843* 

0.795 

278.5 

COMSAT  Non- 
reflective  (CNR) 

0.155 

0.906 

0.803 

0.935 

290.1 

0.155 

0.906 

0.843* 

0.891 

286.6 

*  With  cover  slide  c  =  0.843  (e.g.,  ceria-doped  microcheet). 


IRAS  data  on  the  ANIKs  and  the  WESTARs  show  good  agreement  with  the  model  temperature. 
Satellite  Business  System  Consortium  of  COMSAT,  SBS-1  and  -3  were  reported  to  have  only  10,  instead 
of  20,  20-W  RF  power  TWTs  [24],  With  the  assumption  that  only  about  one-half  the  bus  power  was  used, 
and  the  remaining  power  was  reradiated  at  LWIR  instead  of  RF,  the  CNR  cell  equivalent  of  IJ  was  one- 
half  of  1 5%.  The  result  is  in  good  agreement  with  the  observed  temperature  and  the  solar  cell  model,  thus 
explaining  the  approximate  9°  difference  between  the  295  K  for  the  SBSs  and  the  286  K  for  the  ANIKs 
and  WESTARs. 

The  Galaxy  1  with  K4  %  cells  on  the  fixed  upper  cylinder  and  K7  on  the  deployable  panel  was 
modeled  by  averaging  the  properties  of  the  old  solar  cell  with  the  CNR  cell  properties  from  Table  10.  The 
resulting  model  is  within  about  6  K  of  the  IRAS  measured  temperature. 

The  good  agreement  found  between  the  simple  model  and  the  IRAS  measurements  should  not  lead 
to  the  conclusion  that  the  simple  model  is  an  adequate  treatment  of  a  very  complex  thermodynamic 
system.  Rather,  the  neglected  details,  such  as  the  antenna  thermal  flux  and  radiation  from  the  open 
deployed  cylinder,  tend  to  compensate  for  one  another. 
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TABLE  11 

Physical  Properties  for  HS-376  Cylinders 


SSN 

Satellite 

Model 

Temp 

IRAS 

Temp 

Laboratory 

e 

IRASc 

Laboratory 

a-n 

IRAS 

(a-n)* 

13431 

ANIK  D1 

286.6 

277.1 

0.843 

0.973 

0.751 

0.763 

13652 

ANIK  C3 

286.6 

283.6 

0.843 

0.892 

0.751 

0.768 

13069 

WESTAR  4 

286.6 

281.4 

0.843 

0.891 

0.751 

0.744 

13269 

WESTAR  5 

286.6 

286.4 

0.843 

0.917 

0.751 

0.822 

12065 

SBS-1 

293.7 

296.9 

0.843 

0.728 

0.828 

0.747 

13651 

SBS-3 

293.7 

294.5 

0.843 

0.818 

0.828 

0.815 

14158 

Galaxy  1 

280.1 

278.5 

0.843 

0.836 

0.688 

0.670 

*  (a-n)  is  tabulated  as  a  in  Table  9  IRAS  measurement. 

The  previous  analysis  concentrated  on  IRAS  detections  of  the  eight  HS-376  spin-stabilized 
cylinders.  The  moderate  consistency  of  flux  densities  and  temperatures  encouraged  analysis  of  the 
remaining  observations  of  spin-stabilized  cylinders.  Table  1 2  gives  some  physical  information  about 
many  of  the  satellites  in  orbit  during  IRAS  observations  as  well  as  the  spin-stabilized  cylinders  observed: 
size,  weight,  electrical  power,  type,  name,  and  manufacturer.  Table  1 3  gives  the  observed  flux  densities — 
uncorrected  for  color  temperature — and  the  calculated  temperatures,  emissivities,  and  absorptivities.  All 
other  observed  cylinders  had  smaller  projected  areas  than  the  HS-376  class,  and  almost  all  the  reported 
flux  densities  were  smaller  than  those  for  the  HS-376s.  Two  INTELSAT  4s  (4881  and  6052)  and  three 
INTELSAT  4As  (8620,  10557,  and  10778)  were  observed  to  have  temperatures  between  291  and  298  K, 
with  a  mean  near  295  K.  Both  classes,  HS-312  and  HS-381,  have  approximately  the  same  projected  area. 
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TABLE  12 

Payloads  Detected  by  IRAS 
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TABLE  12  (Continued) 
Payloads  Detected  by  IRAS 
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TABLE  13 

IRAS  Detections  of  Spin-Stabilized  Cylinders 
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The  three  INTELSAT  4A  satellites  with  nearly  equal  flux  densities,  absorptivities,  and  temperatures 
were  selected  as  a  baseline  for  estimating  the  projected  area  and  solar  power  from  the  12-//  flux  density. 
The  mean  value  for  •  -777  establishes  A^-  6.74/1 .777  rr^/J Using  projected  solar  cell  area,  = 
6.74  m^,  the  mean  e  is  computed  to  be  0.928,  and  the  mean  a  is  0.937.  Thus,  the  ale  is  very  close  to 
unity,  and  the  observed  temperature,  295.9  K,  is  very  close  to  that  of  the  free-space  cylinder  model,  which 
is  295  K.  The  INTELSAT  4A  series  provides  a  better  baseline  than  the  newer  HS-376  because  the 
analysis  is  simplified:  the  INTELSAT  series  did  not  use  the  new  K7  high-efficiency  solar  cells.  Table  14 
gives  the  COMSAT  model  temperature  for  the  INTELSAT  4  solar  panels. 


TABLE  14 

INTELSAT  4A  Model  Temperature 


Summer 

Solstice 

Winter 

Solstice 

Equinox 

Eclipse 

Solar  Power 

FWD 

287  K 

290  K 

295  K 

220  K 

AFT 

283  K 

290  K 

293  K 

194  K 

Sun  Shield 

DISC 

285  K 

255  K 

269  K 

233  K 

CONE 

264  K 

244  K 

254  K 

222  K 

TWT  Power 
Supply 

309  K 

312  K 

308  K 

283  K 

310  K 

314  K 

309  K 

276  K 

Figure  1 4  plots  the  measured  (uncolor  corrected)  12-//  flux  vs.  for  the  cylinders  observed  (from 
the  data  in  Table  1 3).  Also  shown  are  parameter  lines  for  1 .0  and  0.5,  which  are  valid  for  temperatures 
near  295  K.  Note  that  most  of  the  satellites  fall  between  e  =  1.0  and  0.75,  although  several  satellites, 
smaller  than  3  m",  give  flux  densities  above  the  €=■  1  parameter  line.  Since  €  cannot  exceed  unity,  the 
implication  is  that  additional  satellite  surface  area  is  contributing  observable  flux.  This  could  be  the 
antenna  or  a  top-of-the-cy Under  optical  solar  radiator.  Equally  likely  is  that  the  sum  of  the  small  signal 
and  the  space  background  yields  excessive  flux  observations.  The  fact  that  many  more  objects  are  above 
the  unity  parameter  line  than  below  suggests  that  the  SNR  alone  is  not  the  cause.  The  detection  threshold 
level,  0.15  Jy,  is  also  shown  on  Figure  14.  One  of  the  smallest  satellites,  Ayame  2,  about  1.2  m^,  has  a 
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12-/X  flux  slightly  greater  than  the  threshold.  Most  of  the  satellites  smaller  than  a  few  square  meters 
yielded  only  one  of  two  detections  compared  to  more  than  a  dozen  for  the  largest.  Some  large  cylinders 
also  have  few  or  a  single  detection  such  as  COMSTAR  2,  =  6.9  m^,  which  appears  to  yield  inaccurate 

results,  e.g.,  T  =  8000  K! 


AREA  (m2) 


Figure  14.  Flux  density  vs.  projected  area. 


As  might  be  expected,  the  temperature  measurement  suffers  greatly  with  small  flux  densities  and 
small  samples.  Recall  the  formal  flux  density  uncertainty  in  the  12-/1 HS-376  flux  data  is  0.6  Jy.  Figure  15 
plots  calculated  temperature  vs.  the  projected  solar  cell  area  for  the  cylinders  observed.  Almost  all  of  the 
largest  satellites  are  at  or  below  300  K.  For  smaller  cylinders,  the  temperature  deviation  from  300  K 
increases.  Arrows  are  shown  for  10%,  20%,  30%,  and  50%  change  in  the  J12/J25  flux  ratio.  Above 
several  square  meters,  most  satellites  are  close  to  or  somewhat  below  295  K. 

The  projected  area  of  a  geostationary  spin-stabilized  satellite,  as  well  as  the  solar  power,  can  be 
estimated  directly  from  the  flux  density.  The  12-/X  flux  (uncolor  corrected)  is  multiplied  by  a  constant  for 
area  and  a  constant  for  power.  The  values  of  the  constants  for  satellites  of  the  INTELSAT  4,  4A  epoch 
and  earlier  are  given  in  Table  15.  Values  for  power  estimation  are  shown  for  beginning  of  life  (BOL)  as 
well  as  end  of  life  (EOL).  The  solar  power  is  based  on  520  W  available  at  EOL,  which  is  seven  years. 
Values  are  also  given  for  the  HS-376  class  and  for  the  somewhat  different  Hughes  Galaxy  1 . 
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249378-15 


Figure  15.  Temperature  vs.  projected  area. 


TABLE  15 

Estimators  for  Geosynchronous  Area  and  Power 


Older  Cells 

HS-376 

Galaxy  1 

INTELSAT  4,4A 

K7  Cells 

K4%  &  K7  Cells 

Area/Jq’^*' 

3.79  m^/Jy 

4.65  m^Jy 

5.48  m^/Jy 

BOL  Power/Jq'2‘' 

365  W/Jy 

512  W/Jy 

614  W/Jy 

EOL  Power/Jq’2‘' 

292  W/Jy 

41 0  W/Jy 

491  W/Jy 

Table  1 6  applies  the  simple  flux  density  relationship  to  a  number  of  satellites  including  Lincoln 
LES-6  (3431),  launched  in  1968.  The  area  and  power  estimates  for  the  GOES  satellites  L  2,  3,  and  4  are 
also  given  together  with  the  Japanese  Himawari  1  and  2.  (All  are  components  of  the  global  Earthwatch 
weather  system  that  forecasters  display  on  television.)  GOES  1  and  synchronous  meteorological  satellite 
(SMS-3)  are  the  series  built  by  Philco  Ford  for  NASA  before  Hughes  began  the  GOES  series.  The  area 


47 


calculated  for  Table  16  agrees  with  the  drawings:  the  power  estimate  has  not  been  verified.  Both 
COMSAT  and  MARISAT  appear  to  yield  good  power  and  area  agreement.  COMSTAR  may  have  some 
new  solar  cells  that  yield  additional  power. 

Most  of  the  smaller  cylinders  are  older  than  the  INTELSAT  4  series  and  should  be  near  the  SOO-K 
temperature.  If  there  is  no  automatic  shutdown,  power  on  the  order  calculated  may  still  be  available  since 
the  degradation  rate  is  slowed  after  ten  years  (5%  in  the  first  year,  18%  after  ten  years). 


TABLE  16 

Geosynchronous  Area  and  Power  Estimation 


Projected  Area  (m^) 

Power  (W) 

Sateiiite 

J 

*'y 

IRAS 

Drawing 

IRAS 

Raferancad 

LES-6  3431 

0.375 

1.43 

1.72 

110 

220  BOL 

190  1  Yr 

GOES  1  8366 

0.756 

2.87 

2.56 

221 

* 

GOES  2  10061 

0.742 

2.82 

2.54 

216 

* 

GOES  4  11964 

0.714 

2.71 

3.35 

209 

* 

GOES  6  14050 

0.732 

2.78 

3.35 

214 

* 

Himawril  10143 

0.738 

2.80 

3.35 

214 

225 

Himawri2  12677 

0.616 

2.34 

2.35 

180 

225 

COMSTAR  4  12309 

1.97 

7.47 

6.9 

575 

760  BOL 

608  7Yr 

Mahsat  3  8882 

0.882 

3.35 

3.45 

258 

220  BOL 

264  EOL 

Marisat2  9478 

0.827 

3.14 

3.45 

242 

220  BOL 

264  EOL 

*  Not  found 

48 


Two  important  thresholds  exist  for  space  surveillance.  The  first,  and  lower,  is  the  minimum  signal 
required  for  detection  for  a  metric  position  measurement.  The  second,  requiring  greater  signal  strength,  is 
the  threshold  for  useful,  repeatable  space  object  information  such  as  emissivity  area  product  and 
temperature.  The  data  in  Table  13  shows  that  cylinders  smaller  than  about  3  m^  usually  exhibit  flux 
densities  below  0.6  Jy,  the  formal  uncertainty  determined  in  the  HS-376  analysis.  Low  flux  densities 
correspond  to  low  SNRs,  hence  derived  temperatures  far  from  300  K  and  emissivities  near  0  or  greater 
than  1,  and  unlikely  values  for  a. 

The  demonstrated  IRAS  geosynchronous  satellite  threshold  is  slightly  below  1  m^.  The  threshold 
for  repeatable  space  object  information  appears  to  be  slightly  above  the  3-m^  cylinder  projected  area 
corresponding  to  about  0.5  Jy.  However,  the  information  threshold  is  strongly  dependent  on  the  number 
of  observations. 
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10.  ROCKET  BODIES 


There  are  668  IRAS  detections  of  rocket  bodies.  The  common  denominator  among  rocket  bodies  is 
their  white  paint  (Ti02)  covering,  which  helps  to  maintain  a  low  temperature  while  on  the  launch  pad.  As 
discussed  in  Section  3,  the  low  temperature  results  from  the  small  value  of  a  /f  =  0.22  for  Ti02.  The  IRAS 
data  on  rocket  bodies  is  represented  in  Figure  16,  which  displays  the  f\2jjf25n  temperature  plotted  against 
the  temperature  for  all  individual  hits.  Figure  17  displays  the  107  detections  on  satellites  for  which 

there  were  more  than  two  detections  (n  >  2),  There  are  three  notable  features.  First  is  the  large  scatter  in 
both  temperatures.  The  scatter  of  these  data  is  significantly  greater  than,  say,  the  spin-stabilized  cylinders. 
A  contributing  factor  could  be  the  slow  passage  across  the  focal  plane  (=  10  sec)  of  a  tumbling  target.  The 
change  in  projected  area  with  time  can  introduce  error  in  the  deduced  flux  ratio.  The  extreme  values  are  as 
low  as  50  K  and  as  high  as  800  K.  Second  is  the  value  of  the  median  temperature  (from  the 
n  >  2  data)  =  182  K.  One  would  expect  (a  /e)*^  x  295  K  =  202  K  for  a  tumbling  cylinder 

covered  with  Ti02.  This  is  consistent  with  the  white  paint  covering.  Third  is  the  difference  between 
T\f\2^f2^^)  =  182  K  and  =  235  K.  The  median  temperatures  for  all  the  detections  are  283  and 

249  K,  respectively.  This  suggests  that  the  emissivity  of  white  paint  has  some  wavelength  dependence. 
This  is  similar  to  the  general  result  (Section  7),  specifically  with  the  results  from  spin-stabilized  cylinders 
(Section  8),  and  the  discussion  in  Section  4 — but  in  the  opposite  sense.  Finally,  Figure  18  gives  the 
temperature  as  a  function  of  SSC  number,  and  Figure  19  gives  the  same  plot  for  satellites  with  more  than 
two  detections.  This  is  the  order  of  launch:  smaller  SSC  numbers  are  associated  with  earlier  launches.  A 
change  in  the  temperature  is  not  apparent.  Dow’s  suggested  increase  in  a  with  time  would  result  in  higher 
temperatures  for  earlier  launches.  This  is  not  observed  in  the  IRAS  data. 
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Figure  16.  Rocket  body  temperatures. 


Figure  17.  Rocket  body  temperatures,  n  >  2. 
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Figure  J 8.  Rocket  body  temperature  vs,  SSC  number. 
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Figure  19.  Rocket  body  temperature  vs.  SSC  number,  n  >  2. 
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11.  THREE-AXIS  STABILIZED  SATELLITES 


11.1  GPS 

The  Global  Positioning  System  (GPS)  provides  a  means  of  determining  position.  These  three-axis 
stabilized  satellites  are  in  12-h  circular  orbits.  The  Block  I  series.  Figure  20  [28],  have  solar  panels 
furnishing  400  W  of  power  and  are  in  orbits  with  an  inclination  of  63.5°.  The  Block  II  series,  Figure  21, 
have  solar  panels  furnishing  700  W  of  power  and  are  in  orbits  with  an  inclination  of  56°. 

The  Block  I  solar  panels  are  in  the  shape  of  a  cylindrical  section,  whereas  the  Block  II  solar  panels 
are  flat.  The  satellite  attitude  is  continuously  adjusted  so  the  antenna-populated  side  of  the  satellite  is 
facing  the  earth,  and  the  solar  panel  arrays  are  normal  to  the  sun  direction.  For  the  Block  I  satellites,  the 
earth-facing  side  has  an  area  of  =  1.35  m^,  and  the  solar  panels  have  a  projected  area  of  = 
5.28  m^.  To  complete  the  model,  assume  the  main  body  of  a  GPS  satellite  is  cylindrical  in  shape. 
Therefore,  the  side  silhouette  is  a  rectangle  of  area  =  1 .6  m^  Further,  assume  that  this  cylinder  and 
flat  solar  panel  model  is  a  simple  model  that  has  no  shadowing. 

The  IRAS  satellite  made  22  detections  of  five  Block  I  satellites:  SSC  numbers  10684,  11054, 
1 1141,  1 1690,  and  1 1783.  These  measurements  were  made  with  sun-GPS-IRAS  phase  angles  ranging 
from  78°  to  120°.  Temperatures  determined  with  the  /12//251J  ^ux  ratio  are  plotted  in  Figure  22  as  a 
function  of  time.  There  are  some  significant  outliers.  The  median  temperature  is  297  K  and  is  preferred  to 
the  average  to  avoid  a  bias  from  the  extreme  values.  Temperatures  determined  with  the  f^^x  ratio 

are  plotted  in  Figure  23  as  a  function  of  phase  angle.  Except  for  three  outliers,  the  temperature  seems  to 
be  independent  of  phase  angle.  This  suggests  that  the  solar  panels  have  the  same  temperature  as  the  body. 
The  reported  1 2-fi  flux  as  a  function  of  phase  angle  is  plotted  in  Figure  24.  A  dependence  of  flux  on  phase 
angle  is  shown.  Averaging  groups  of  observations,  one  observes  a  minimum  of  1.07  Jy  at  d  =  90.6°, 
1.57Jy  at0=  101.8°,  raising  to  2.71  Jy  at  6=  120.8°  and  1.72  Jy  at  0=81.3°. 
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249378-20 


Figure  20.  GPS  Block  I  satellite. 


56 


249378-21 


Figure  21.  GPS  Block  1 1  satellite. 
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Figure  22.  GPS  temperature  V5.  date. 
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Figure  23.  GPS  temperature  vs.  phase  angle. 
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74U7A-24 
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Figure  24.  GPS  l2-jl  flux  V5.  phase  angle. 


As  an  illustration,  the  emissivity  is  calculated  assuming  the  0  =  90®  data  all  from  the  earth-facing, 
1.35  m^  side  of  the  GPS  satellite  is  the  only  source  of  radiant  energy.  So  the  radiant  flux  density  at  300  K 
is 


=  27.39  W/\i/  (effective  eA) 


The  IRAS  color  correction  adjusts  the  90®  value  to 


ji2ix  ^  =  1.07/0.92  =  1.163  Jy 

obs  g  12u 


which  is  equal  to  2.421  x  10^**^  W/(mV).  Now  at  a  typical  range  of  19,000  km,  the  observed 

spectral  radiation  intensity  is  k{\9.0  x  10^)^(2.421  x  10“'*^)  =  27.46  W//i.  Therefore,  the 

effective  area  eA  is 


eA= 


nr2pi2u 

obs 


1.0025  (m^) 
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Since  the  projected  area  is  A  =  1 .35  m^,  the  emissivity  of  the  flat,  antenna-populated  side  is  e  =  0.74. 

A  more  detailed  model  can  be  defined  as  follows.  At  phase  angles  6  ^  90°  and  aspect  angles  <p  ^0, 
the  solar  panel  and  the  satellite  body,  assumed  cylindrical,  will  contribute  to  the  observed  flux.  Now  the 
cylinder  and  plate  model  will  have  the  projection  of  the  solar  panel  on  the  IRAS  line  of  sight  as 
Ajoii]jlcos(0)l.  So  if  the  cylinder  flat  end  area  is  1.35  m^,  the  cylinder  side  projected  area  is  = 

1 .6  m^,  and  the  solar  panel  area  is  A  5.28  m^,  then  the  flux  received  would  be 


P  ° 
^  obs 


llR 


(cp)  J  S  in  ( (p )  |  )  .^,  J  J  COS  (  0  )  | 


where  (j)  is  the  IRAS-GPS  aspect  angle.  The  back  of  the  solar  panel  is  observed  at  0  >  90®  and  the 
emissivity  is  e";  the  front  of  the  solar  panel  is  observed  at  0  <  90®  and  the  emissivity  is  The  illustrative 
calculation  can  be  generalized  to  obtain  the  emissivity  of  the  solar  panels  accounting  for  the  phase  angle, 
the  aspect  angle,  and  the  range.  The  linear  regression  results  are  summarized  in  Table  17. 


TABLE  17 

Emissivity  for  GPS  Sateliites 


0 

Element 

AK) 

e 

=90 

Body  Flat 

1.35 

0.622±0.04 

<90 

Solar  Cell  Front 

5.28 

0.933±0.10 

>90 

Solar  Cell  Back 

5.28 

0.619±0.08 

The  emissivity  for  the  solar  panel  front  side,  e  =  0.93  ±0.l,  is  larger  than  expected.  It  depends 
directly  on  the  adopted  area.  Based  on  laboratory  measurements  one  expects  emissivity  values  about  0.80. 
The  spin-stabilized  cylinders  give  values  of  solar  cell  emissivity  between  0.6  and  0.7.  The  solar  panel 
back  and  the  main  body  have  the  same  emissivity,  e  =  0.62.  The  two  values  are  equal  within  the  statistical 
uncertainty. 

11.2  GORIZONT 

The  Gorizont  satellites,  for  telephone  and  international  television,  were  launched  by  the  Soviet 
Union.  Gorizont  I  (SSN  1 1 158)  has  a  24-h  period  and  was  launched  in  1978  into  an  orbit  with  1 1.3® 


60 


inclination.  By  1993  the  inclination  increased  to  more  than  21°.  The  other  Gorizont  satellites  are 
geosynchronous  satellites  in  low  inclination  orbits. 

The  Gorizont  is  a  three-axis  stabilized  satellite  (see  Figure  25)  [28].  The  main  body  is  about  5  m 
long  and  2  m  in  diameter.  There  are  a  number  of  appendages  that  will  increase  the  effective  area.  Two 
panels  of  solar  cells  of  unknown  size  provide  power. 


Figure  25.  Gorizom  satellite.  (Used  with  permission.  Donald  H.  Martin,  Communication  Satellites  1958-1992, 
©1991,  The  Aerospace  Corporation.) 


The  IRAS  satellite  made  26  detections  of  four  Gorizont  satellites:  SSC  numbers  1 1 158,  1 1440, 
13092,  and  13624.  These  measurements  were  made  with  sun-Gorizont-IRAS  phase  angles  ranging  from 
60°  to  1 17°.  Temperatures  determined  with  the  ratio  are  plotted  in  Figure  26  as  a  function  of 

time.  The  median  temperature  is  280  K  and  is  preferred  to  the  average  to  avoid  a  bias  from  the  extreme 
values.  Temperatures  determined  with  the/, 2^/25^  flux  ratio  are  plotted  in  Figure  27  as  a  function  of  phase 
angle.  The  trend  is  temperature  independent  of  phase  angle.  This  suggests  that  the  solar  panels  have  the 
same  temperature  as  the  body.  The  reported  12-//  flux  as  a  function  of  phase  angle  is  plotted  in  Figure  28. 
Here,  a  minimum  of  1.3  Jy  is  observed  at  0  =  90°,  raising  to  2.2  Jy  at  0  =  80°,  and  2.8  Jy  at 
0  =  115°.  Note  also  that  none  of  the  data  for  Gorizont  I  follow  this  trend.  Therefore,  it  is  not  included  in 
the  following  calculation. 
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Figure  26.  Gorizont  temperature  v^.  date. 
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Figure  27.  Gorizont  temperature  vs.  phase  angle. 
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Figure  28.  Gorizonl  I2-fl  flux  vs.  phase  angle. 


Now  the  radiant  flux  density  at  280  K  is 


^.^^^{280)  =  21.06  W/]ji/  {effective  eA) 


The  IRAS  color  correction  adjusts  the  90°  value  to 

J-12U  ^  _  1.3/0.90  =  1.44  Jy 

ohs  q  12u 


which  is  equal  to  -  3.007  x  lO"'**  W/(m^/i).  Now  at  a  range  of  36,000  km,  the  observed  spectral 

radiation  intensity  is  n^F^'^  =  7r(36.0  x  10*)^(3.007  x  10"''*)  =  122.4  W///.  Now  the  effective  area  eA  is 


eA= 


nr^F 


12u 

obs 


•^12u(280) 


=  5.81 


In  this  case  we  chose  to  adopt  a  value  of  the  emissivity,  e  =  0.70,  and  calculate  the  earth-facing  area, 
which  is  A  =  8.3  m^.  This  is  equivalent  to  a  3.2-m-diameter  main  body,  which  is  60%  larger  than  the 
nominal  main  body  dimension.  We  can  estimate  the  size  of  the  solar  panels.  Following  the  same 
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reasoning  used  for  GPS  (see  Section  1 1.1),  the  projection  of  the  solar  panel  on  the  IRAS  line  of  sight  is 
Ajoiaflcos(0)l.  So  if  the  cylinder  flat  end  area  is  /!(,»,,  the  cylinder  side  project  area  is  and  the  solar 
panel  area  is  A^,  then  the  flux  received  would  be 

((P)  +A^,^Jsin((p)  I  (0)  ] , 

where  (p  is  the  IRAS-Gorizont  aspect  angle.  The  back  of  the  solar  panel  is  observed  at  0=115°,  and  the 
front  is  observed  at  0  =  80°.  The  same  regression  calculation  can  be  done.  Since  we  have  no  prior 
knowledge  of  the  solar  plane  area,  we  chose  to  determine  the  emissivity  area  product  eA.  The  results  are 
summarized  in  Table  18. 


TABLE  18 

Emissivity  Area  for  Gorizont  Sateliites 


e 

Element 

eA  (m*) 

A(m^) 

=90 

Body  Flat 

4.38±0.6 

6.25 

<90 

Solar  Cell  Front 

30.75±6.8 

33 

>90 

Solar  Cell  Back 

15.50±2.2 

25 

This  estimate  of  the  earth-facing  area,  eA  =  4.38,  leads  to  an  area  A  =  6.25  m^,  equivalent  to  a  main 
body  diameter  of  2.8  m,  which  is  about  40%  larger  than  the  nominal.  The  estimates  of  the  panel  size  from 
the  front  and  back  emissions  differ  due  to  measurement  error,  model  error,  and  above  all,  error  in  the 
adopted  emissivity.  Based  on  the  GPS  results  one  might  expect  the  back  panel  emissivity  to  be  lower  by 
30%  to  50%.  If  we  adopt  the  GPS  values  of  emissivity  from  the  body  and  solar  cells,  then  we  get  values 
given  in  Table  1 8.  The  solar  panel  area  is  of  the  order  of  25  to  35  m^.  This  is  consistent  with  two  solar 
panels,  each  providing  1  to  2  kW  of  power. 

11.3  RADUGA 

The  Raduga  satellites  were  the  first  Soviet  satellites  in  geosynchronous  orbit.  They  were  three-axis 
stabilized  and  were  used  for  telephone,  telegraph,  and  television  transmission.  The  characteristics  are 
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similar  to  the  Gorizont  satellites  (Section  1 1.2).  The  IRAS  satellite  made  27  detections  of  six  Raduga 
satellites:  SSC  numbers  10159,  10987,  1 1708,  12897,  13669,  and  13974.  These  measurements  were  made 
with  the  sun-Raduga-IRAS  phase  angle  ranging  from  84°  to  103°.  Temperatures  determined  with  the 
/i2^/25P  plotted  in  Figure  29  as  a  function  of  time.  These  measurements  were  made  with  sun- 

Raduga-IRAS  phase  angles  ranging  from  76°  to  103°.  Temperatures  determined  with  the /i 2^/25^  flux  ratio 
are  plotted  in  Figure  30  as  a  function  of  phase  angle.  The  median  temperature  is  288  K  and  is  preferred  to 
the  average  to  avoid  a  bias  from  the  extreme  values.  The  trend  is  temperature  independent  of  phase  angle. 
This  suggests  that  the  solar  panels  have  the  same  temperature  as  the  body.  The  reported  12-/Z  flux  as  a 
function  of  phase  angle  is  plotted  in  Figure  31.  The  same  regression  calculation  that  was  done  for  GPS 
and  Gorizont  satellite  data  can  be  done.  Since  there  is  no  prior  knowledge  of  the  solar  plane  area,  we 
determined  the  emissivity  area  product  eA.  The  results  are  summarized  in  Table  19. 
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Figure  29.  Raduga  temperature  vs.  date. 
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Figure  30,  Raduga  temperature  v^.  phase  angle. 
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Figure  3L  Raduga  I2-pflux  vs,  phase  angle. 
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TABLE  19 

Emissivity  Area  for  Raduga  Satellites 


0 

Element 

€A  (m^) 

A(m^) 

=90 

Body  Flat 

2.80±1.4 

4.00 

<90 

Solar  Cell  Front 

37.13±14.4 

40 

>90 

Solar  Cell  Back 

36.44±16.4 

58 

This  estimate  of  the  earth-facing  area,  eA  =  2.8  m^,  leads  to  an  area  >4  =  4  m^,  which  is  equivalent  to 
a  main  body  diameter  of  2.25  m.  There  is  no  prior  knowledge  of  the  body  size.  The  estimates  of  the  panel 
size  from  the  front  and  back  emissions  differ  due  to  measurement  error,  model  error,  and  above  all,  error 
in  the  adopted  emissivity.  Based  on  the  GPS  results  one  might  expect  the  back  panel  emissivity  to  be 
lower  by  30%  to  50%.  However,  the  e/\  product  for  the  front  and  back  are  nearly  equal.  Nevertheless, 
uncertainty  is  quite  large.  If  the  GPS  values  of  emissivity  from  the  body  and  solar  cells  are  adopted,  then 
values  for  A  are  as  given  in  Table  1 8.  The  solar  panel  area  is  larger  than  35  m^.  This  is  consistent  with  two 
solar  panels,  each  providing  1  to  2  kW  of  power. 

11.4  EKRAN 

The  Ekran  satellites  were  three-axis  stabilized  and  were  used  for  telephone,  telegraph,  and 
television  transmission  within  the  USSR  and  bordering  countries.  The  characteristics  are  similar  to  the 
Gorizont  satellites  (Section  1 1.2),  Figure  32  [28],  with  the  addition  of  a  complex  antenna  array  on  the 
earth-facing  side  of  the  satellite.  The  IRAS  satellite  made  27  detections  of  five  Ekran  satellites:  SSC 
numbers  1 1273,  12120,  12564,  13554,  and  13878.  These  measurements  were  made  with  the  sun-Ekran- 
IRAS  phase  angle  ranging  from  70°  to  97°.  Temperatures  determined  with  the  flux  ratio  are 

plotted  in  Figure  33  as  a  function  of  time.  Temperatures  determined  with  the  fxijf25^  ^ux  ratio  are  plotted 
in  Figure  34  as  a  function  of  phase  angle.  The  median  temperature  is  284  K  and  is  preferred  to  the  average 
to  avoid  a  bias  from  the  extreme  values.  There  is  a  suggestion  of  an  increase  in  temperature  for  phase 
angles  greater  than  90°.  The  reported  12-/i  flux  as  a  function  of  phase  angle  is  plotted  in  Figure  35.  The 
temperature  measurements  are  consistent  with  those  of  other  satellites.  However,  there  seems  to  be  no 
simple  relation  between  phase  angle  and  flux.  A  more  complex  radiation  model  is  necessary  to  explain 
these  data,  which  are  probably  due  to  the  large  antenna  array. 
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249378-38 


Figure  32.  Ekran  satellite.  (Used  with  permission.  Donald  H.  Martin.  Communication  Satellites  1958—1992, 
©1991,  The  Aerospace  Corporation.) 


11.5  ATS-6 

The  Application  Technology  Satellite  (ATS)-6,  SSC  number  7318,  was  active  for  more  than  five 
years.  It  was  used  for  a  number  of  experiments  in  satellite  communications  technology.  It  was  a  three-axis 
stabilized  satellite  with  an  enormous  parabolic  antenna,  a  split  cylindrical  solar  cell  panel,  and  a 
rectangular  payload;  see  Figure  36  [28]. 

The  IRAS  satellite  made  three  detections  of  the  ATS-6,  which  are  given  in  Table  20.  The  curiosity 
is  the  large  flux  measured.  The  antenna  material  was  copper  netting  on  a  structure  of  48  aluminum  ribs,  of 
which  there  is  no  specific  knowledge  of  the  emissivity  or  temperature.  Still,  even  with  extreme  values  of 
the  emissivity,  it  is  difficult  to  calculate  such  a  large  flux. 
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Figure  33.  Ekran  temperature  vs.  date. 


PHASE  ANGLE  (deg) 

Figure  34.  Ekran  temperature  v.y.  phase  angle. 
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^2-\l  FLUX  (Janskys) 


Figure  35.  Ekran  12-/J.  flux  vs.  phase  angle. 
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24937»^2 


Figure  36.  ATS-6  satellite.  (Used  with  permission.  Donald  H.  Martin,  Communication  Satellites  1958-1992, 
©1991,  The  Aerospace  Corporation,) 
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TABLE  20 

ATS-6  Observed  Flux 


MJD 

J 

''a 

'ri2ii/25u 

T25u«Ou 

0 

4) 

45521.605 

4.88 

6.98 

2.24 

284 

349 

92 

0.3 

45551.092 

5.48 

7.45 

2.24 

291 

379 

96 

1.8 

45641.552 

5.51 

9.00 

3.42 

267 

287 

89 

0.8 

72 


12.  DISCUSSION 


The  IRAS  satellite  made  both  position  and  radiometry  observations  of  RSOs.  The  accuracy  of  the 
position  measurement  was  limited  by  the  detector  size  of  4.5  arcmin.  This  was  adequate  for  correlation 
with  a  satellite  catalogue  to  associate  observations  with  specific  satellites.  Use  of  these  position 
observations  for  determination  of  ephemerides  is  not  addressed  here.  Concerning  the  radiometry 
measurements,  the  accuracy  of  a  single  observation  is  experimentally  determined  to  be  ±0.6  Jy.  Of  the 
190,000  detections,  2047  were  correlated  with  the  known  satellite  catalogue.  The  correlation  used  both 
the  position  and  velocity  measurement.  The  observations  were  correlated  with  452  satellites.  Due  to  the 
zenith  pointing  of  IRAS,  most  of  the  observations  are  on  deep  space  satellites,  many  of  which  are 
geosynchronous. 

The  radiometer  flux  data  required  some  screening.  The  principal  issue  concerned  partial  detections 
where  the  image  crossed  the  detector  edge.  Lacking  detailed  knowledge  of  the  image  path  across  the 
detector,  this  could  only  be  done  statistically.  Therefore,  even  the  screened  data  will  have  some  partial 
detections.  Nevertheless,  screening  individual  detections  significantly  improved  the  data  consistency. 

Calibration  of  the  IRAS  data  remains  a  nagging  question.  There  have  been  suggestions  about 
revision  of  the  IRAS  calibration,  as  reviewed  in  Section  5.  There  are  significant  questions  of 
interpretation  that  would  be  affected  by  a  change  in  calibration.  For  example,  for  all  the  IRAS  data  there 
is  the  uniformly  lower  temperature  observed  with  the  25-  to  60-//  band  flux  ratio  than  observed  with  the 
12-  to  25-//  band  flux  ratio.  This  is  true  for  the  majority  of  the  observations  and  is  particularly  striking  in 
the  Hughes  HS-376  cylinders.  As  Table  3  shows,  most  critiques  suggest  that  the  25-//  fluxes  are  high  by 
4%  to  10%.  For  our  analysis  of  the  HS-376  cylinders,  changing  the  calibration  in  this  way  would  result  in 
worse  agreement  between  the  IRAS  observations  and  the  models.  Also,  any  decrease  in  the  25-//  flux 
would  increase  the  disagreement  between  the  12-///25-//  and  the  25-///60-//  temperatures  and  would 
increase  the  derived  emissivities.  On  the  other  hand,  for  both  classes,  the  HS-376  and  the  INTELSAT  4 
and  4A  satellites,  the  average  value  measured  for  emissivity  is  Just  below  unity.  This  is  5%  to  15%  above 
the  laboratory  measurements  for  solar  cells.  Perhaps  this  is  due  in  part  to  the  simple  model  ignoring  the 
antenna  heat.  Perhaps  some  of  the  excess  flux  density  should  be  attributed  to  calibration  error  on  the  order 
of  5%  too  high  in  both  the  12-  and  25-//  bands.  How  to  apportion  these  disagreements  between 
dependence  of  solar  cell  emissivity  on  wavelength  and  calibration  remains  undecided.  Therefore,  no 
definite  conclusion  about  IRAS  calibration  can  be  reached  from  this  analysis. 

IRAS  observations  were  serendipitous.  Therefore,  there  is  no  systematic  set  of  observations  for 
analysis.  On  average  there  were  six  observations  for  each  satellite,  but  many  satellites  were  observed  only 
once.  This  was  ameliorated  by  using  classes  of  satellites  that  can  be  assumed  to  be  virtually  identical,  for 
example,  the  Hughes  HS-376  spin-stabilized  cylinders,  the  GPS  satellites,  and  a  number  of  Soviet 
communication  satellites.  In  this  case  a  data  set  exceeding  20  detections  could  be  formed. 

The  discussion  of  IR  satellite  radiometry  inevitably  turns  on  the  properties  of  spacecraft  materials. 
These  include  paints,  solar  cells,  antennas,  and  radiators.  Existing  laboratory  data  is  suggestive  but  not 
conclusive.  Also,  data  on  the  effects  of  the  space  radiation  environment  on  spacecraft  materials  from  the 
LDEF  are  helpful.  Nevertheless,  the  IRAS  data  must  stand  on  its  own.  For  example,  the  effects  of  Ti02  on 
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rocket  body  temperature  is  consistent  with  expectations.  However,  expectations  that  TiOj  will  change 
absorptivity  after  exposure  to  the  space  environment  are  not  supported  by  the  IRAS  data. 

With  these  caveats,  one  can  see  a  lot  of  information  in  the  IRAS  radiometry.  Even  with  very  simple 
models  one  can  learn  about  the  absorptivity  and  emissivity  of  spacecraft  materials.  Alternatively,  by  using 
the  geometrical  variation  of  phase  angle  and  aspect  angle,  one  can  estimate  the  physical  size  of  certain 
spacecraft  components.  Further  progress  will  entail  considerably  more  complex  models  and  detailed 
knowledge  of  the  spacecraft. 

The  fact  that  the  IRAS  was  a  space-based  platform  was  extremely  valuable.  The  whole  question  of 
atmospheric  absorption  is  avoided,  and  the  data  have  remarkable  consistency.  In  addition,  the 
multispectral  data  immediately  provides  information  on  the  object  temperature  and  leads  directly  to 
physical  analysis.  Much  analysis  of  the  ground-based  data  has  suffered  from  both  these  limitations. 

In  1992  no  data  on  geosynchronous  objects  were  published  from  either  ground-  or  space-based 
LWIR  sensors.  In  the  winter  of  1992-1993,  Seniw  and  Rieke  collected  LWIR  data.  This  included 
60  tracks  on  20  geosynchronous  satellites  in  the  N-  and  Q-bands,  which  allowed  temperature 
measurements  on  two  objects  [3].  The  data  were  obtained  over  a  wide  range  of  sun-sensor-satellite  phase 
angles  and  included  observations  of  temperature  in  and  out  of  the  earth  shadow.  Six  geosynchronous 
objects  observed  by  Seniw  and  Rieke  were  observed  by  IRAS  a  decade  earlier.  Two  cylinders,  GOES  6 
and  the  Hughes  Galaxy  1 ,  had  illumination  phase  angles  of  52“  and  64“,  respectively.  This  is  close  enough 
to  the  IRAS  near  90“  phase  angle  to  warrant  comparison.  Only  astronomical  N-band  data  were  obtained 
from  the  ground  so  that  temperature  could  not  be  determined.  The  radiant  intensity,  watts/steradian,  was 
computed  from  the  IRAS  temperature  and  satellite  eA^  and  multiplied  by  the  fraction  of  the  total  in 
N-band.  Table  21  compares  the  University  of  Arizona  and  the  IRAS  measurements.  In  this  calculation, 
following  the  logic  in  Section  9,  the  EOL  temperature  estimate  for  Galaxy  1  of  295  K  is  used  instead  of 
the  observed  value  of  278  K. 

Recall  that  the  IRAS  pointing  logic  generally  viewed  sun-tracking  solar  panels  nearly  edge-on. 
Therefore,  for  those  satellites,  the  radiant  intensity  was  probably  from  the  equipment  body  and  antennas, 
which  is  a  small  fraction  of  the  RSO’s  radiant  intensity  at  maximum  phase.  Slight  changes  in  aspect  angle 
would  yield  moderate  changes  in  observed  intensity,  which  we  have  exploited.  On  the  other  hand,  spin- 
stabilized  cylindrical  solar  panels  present  a  near-constant  geometry.  The  despun  antennas  are  pointing  at 
the  earth  but  are  edge-on  to  the  sun. 
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TABLE  21 

Ground-Based  and  IRAS  Observations 


Category 

Univ  Arizona 

IRAS 

GOES-6 

Weather  Sat 

14050 

Phase  Angle 

52° 

90o 

Reported  Flux 

0.74  Jy  @  10.6  p 

0.732  Jy  @  12  p 

Temperature 

- 

325  K 

Radiant  Inten. 

-  169  W/sr 

132  W/sr 

Galaxy  1 

HS-376 

14158 

Phase  Angle 

64o 

90o 

Reported  Flux 

1.58  Jy  @  10.6  M 

1.61  Jy  @  12  p 

Temperature 

295  K 

Radiant  Inten. 

372 

382 
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13.  CONUNDRUMS 


The  IRAS  debris  data  base  presents  a  number  of  unanswered  questions — some  going  to  the  heart  of 
this  analysis.  In  Section  9,  analysis  of  Hughes  HS-376  spin-stabilized  cylinders  (the  simplest, 
homogeneous,  and  unambiguous  data  set  available)  is  presented.  The  radiometric  model  is  quite  simple 
and  reasonably  good.  The  data  set,  after  screening,  seems  consistent,  and  the  derived  results  are  quite 
reasonable.  From  the  analysis  of  variance,  we  derived  the  formal  uncertainty  of  a  single  observation  as 
±0.6  Jy  in  the  12-/i  band.  However,  to  reach  this  state,  some  data  had  to  be  edited.  There  were  far  more 
screened  detections,  with  errors  exceeding  1.2  Jy,  than  expected  on  the  basis  of  Gaussian  statistics.  One 
could  postulate  that  another  satellite  was  observed  at  that  time  and  direction.  If  true,  then  where  is  the 
measurement  of  the  correlated  satellite  that  we  know  was-there  at  that  time?  One  must  assume  that  the 
detection  was  indeed  of  the  correlated  satellite,  with  quite  a  large  flux  error.  Consequently,  we  have 
considerable  concern  about  the  error  statistics  of  the  IRAS  debris  data  base.  Unfortunately,  there  is 
insufficient  data  to  actually  discuss  the  error  distribution  function. 

A  data  set  with  a  related  problem  is  derived  as  follows.  Select  all  the  detections  that  were  obtained 
with  small  velocity  (<  6  arcmin/sec)  across  the  focal  plane.  There  are  22,000  such  detections.  Of  these, 
2047  were  correlated  with  the  catalogue.  Excluding  the  pathological  case  of  a  co-orbiting  satellite,  these 
would  be  observations  of  satellites  at  extreme  range:  geosynchronous  satellites  or  high  eccentricity 
satellites  observed  at  apogee.  In  both  cases  the  observed  range  would  be  36,000  km  or  greater.  Among 
these  data,  there  are  54  detections  with  a  12-/Z  flux  >  100  Jy  and  10  exceeding  900  Jy.  Note  that  a  flux 
greater  than  999.99  Jy  is  listed  as  999.99.  Of  these,  only  nine  had  a  12-  to  25-//  flux  ratio  that  implied  a 
plausible  temperature,  e.g.,  between  250  K  and  400  K.  Most  of  the  large  fluxes  occur  in  only  one  band. 
There  are  cases  where  both  detectors  for  a  band  registered  these  high  values.  Only  three  were  observed  at 
small  declination,  i.e.,  in  the  geosynchronous  belt,  with  most  having  very  high  declination.  These  very 
large  flux  detections  seem  to  have  no  consistent  orbits.  One  must  conclude  that  these  are  all  bogus 
observations. 

The  reliability  of  a  single  detection,  even  with  multiple  hits,  is  seriously  in  question.  With  the  spin- 
stabilized  cylinders,  there  were  hits  with  twice  the  “correct”  flux  value.  These  can  be  deduced  because  in 
this  case  there  are  multiple  detections  of  a  well-understood  target.  Detections  at  a  geosynchronous  range 
of  100  Jy  or  more  seem  implausible.  It  is  difficult  to  assess  the  observed  flux  without  multiple  detections. 
Therefore,  we  believe  that  any  conclusions  based  on  a  single  detection  must  be  made  very  cautiously. 
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APPENDIX  A 

IRAS  DETECTIONS  CORRELATED  WITH  THE  KNOWN  SATELLITE 

CATALOGUE 


This  api>endix  is  a  list  of  the  2047  IRAS  detections  correlated  with  the  known  satellite  catalogue. 
This  selection  includes  the  detections  in  agreement  within  0.6°  In  some  cases  detections  were  not  used  in 
further  analysis.  The  reasons  include  angular  velocity  error  exceeding  some  value,  the  radiometric  data 
being  noisy,  or  exceeding  a  statistical  test,  as  discussed  in  the  text. 

The  list  contains  12  items  for  each  detection.  First  is  the  satellite  identification  number.  We  have 
chosen  the  identification  given  by  the  US  Air  Force  Space  Surveillance  Center  (SSC).  In  Appendix  B  we 
provide  other  information  such  as  the  international  designator  and  orbital  characteristics.  All  12  items,  in 
the  following  order,  are  given  below. 

1.  Satellite  identification  number. 

2.  Modified  Julian  date  of  the  detection. 

3.  Temperature,  Kelvin,  derived  from  the  flux  ratio. 

4.  Temperature,  Kelvin,  derived  from  the  flux  ratio.  If  the  flux  ratio  is  not  available,  then  the 

temperature  is  reported  as  -1 . 

5.  Quoted  12- flux,  Janskys.  Data  screening  is  applied  as  described,  eliminating  partial  detections,  and 

the  selected  detections  are  averaged  to  obtain  the  value  given  here. 

6.  Quoted  25- fi  flux,  Janskys,  screened  in  the  same  way  as  the  12-/idata. 

7.  Quoted  60- fi  flux,  Janskys,  screened  in  the  same  way  as  the  \2-fi  data. 

8.  IRAS-satellite-sun  phase  angle,  degrees. 

9.  IRAS  satellite  aspect  angle,  degrees. 

10.  Satellite  declination  of  the  detection. 

1 1 .  IRAS-to-satellite  range  derived  from  the  ephemeris,  km. 

12.  Observed  angular  rate  of  the  detection,  arcmin/sec. 
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APPENDIX  B 

SATELLITES  OBSERVED  BY  IRAS 
(AT  LEAST  ONCE) 


This  appendix  identifies  the  465  satellites  observed  by  IRAS.  This  list  contains  the  following 
information. 

1 .  Air  Force  Space  Surveillance  Center  identification  number. 

2.  Year  of  launch. 

3.  International  designator. 

4.  Country  of  origin. 

5.  Type  of  object:  PL  denotes  a  payload;  RB,  a  rocket  body;  and  DB  is  debris. 

6.  Status  of  the  object. 

7.  Common  name  of  object  or  a  comment. 

8.  Inclination  of  the  object,  degrees. 

9.  Eccentricity  of  the  object. 

10.  Semimajor  axis  of  the  object,  earth  radii. 

1 1 .  Mean  motion  of  the  object,  revolutions  per  day. 
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Sid  Year  Designation  Owner  Mission  Status  Comment  Inclination  Eccentricity  Semi  major  axis  Mean  motion 
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4478  1970  OSSA  ITS  PL  NA  INTELSAT  riTF-8 


Sid  Year  Designation  Owner  Mission  Status  Comment  Inclination  Eccentricity  Semi  major  axis  Mean  motion 
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7170  1973  097D  SOV  RB  NA  MOLNIYA  1-26//S1 


Sid  Year  Designation  Owner  Mission  Status  Coiwnent  Inclination  Eccentricity  Senii  major  axis  Mean  motion 
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Sid  Year  Designation  Owner  Mission  Status  Comment  Inclination  Eccentricity  Semi  major  axis  Mean  motion 
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